Effective action approach to higher-order relativistic tidal interactions in binary 
systems and their effective one body description 
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The gravitational-wave signal from inspiralling neutron-star-neutron-star (or black-hole-neutron- 
star) binaries will be influenced by tidal coupling in the system. An important science goal in 
the gravitational-wave detection of these systems is to obtain information about the equation of 
state of neutron star matter via the measurement of the tidal polarizability parameters of neutron 
stars. To extract this piece of information will require to have accurate analytical descriptions 
of both the motion and the radiation of tidally interacting binaries. We improve the analytical 
description of the late inspiral dynamics by computing the next-to-next-to-leading order relativistic 
correction to the tidal interaction energy. Our calculation is based on an effective- act ion approach 
to tidal interactions, and on its transcription within the effective-one-body formalism. We find that 
second-order relativistic effects (quadratic in the relativistic gravitational potential u — G{mi + 
m2)/{c?r)) significantly increase the effective tidal polarizability of neutron stars by a distance- 
dependent amplification factor of the form 1 + a\u + ot2U^ + ■ ■ ■ where, say for an equal-mass 
binary, ai = 5/4 = 1.25 (as previously known) and 02 = 85/14 ~ 6.07143 (as determined here for 
the first time). We argue that higher-order relativistic effects will lead to further amplification, and 
we suggest a Pade-type way of resumming them. We recommend to test our results by comparing 
resolution-extrapolated numerical simulations of inspiralling-binary neutron stars to their effective 
one body description. 
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I. INTRODUCTION 

Inspiralling binary neutron stars are among the most 
promising sources for the advanced versions of the cur- 
rently operating ground-based gravitational-wave (GW) 
detectors LIGO/Virgo/GEO. These detectors will be 
maximally sensitive to the inspiral part of the GW sig- 
nal, which will be influenced by tidal interaction between 
the two neutron stars. An important science goal in 
the detection of these systems (and of the related mixed 
black-hole-neutron-star systems) is to obtain information 
about the equation of state of neutron-star matter via 
the measurement of the tidal polarizability parameters 
of neutron stars. The analytical description of tidally in- 
teracting compact-binary systems (made of two neutron 
stars or one black hole and one neutron star) has been 
initiated quite recently [ll-ll|- In addition, these analyti- 
cal descriptions have been compared to accurate numeri- 
cal simulations 0, [^-[Hl , and have been used to estimate 
the sensitivity of GW signals to the tidal polarizability 
parameters pll - [l5| . 

Here, we shall focus on one aspect of the analytical 
description of tidally interacting relativistic binary sys- 
tems, namely the role of the higher-order post-Newtonian 



(FN) corrections in the tidal interaction energy, as de- 
scribed, in particular, within the effective one body 
(BOB) formalism Indeed, the analysis of Ref. [Ef, 

which compared the prediction of the EOB formalism 
for the binding energy of tidally interacting neutron 
stars to (nonconformally flat) numerical simulations of 
quasi-equilibrium circular sequences of binary neutron 
stars [20, HH , suggested the importance of higher-order 
FN corrections to tidal effects, beyond the first post- 
Newtonian (IFN) level, and their tendency to signifi- 
cantly increase the "effective tidal polarizability" of neu- 
tron stars. 

In the EOB formalism, the gravitational binding of a 
binary system is essentially described by a certain "radial 
potential" A{r). In the tidal generalization of the EOB 
formalism proposed in Ref. Q , the EOB radial potential 
A{r) is written as the sum of three contributions. 



A{r) = A^^^{r) + AX^\r) + A'i^^'(r) 



tidal / 



1 tidal / 



(1.1) 



where ^^bh ^^le radial potential describing the dy- 
namics of binary black holes, and where A^'^^'(r) and 
A^^^^[r) are the additional radial potentials associated, 
respectively, with the tidal deformations of body A and 
body B. [For binary neutron-star systems both A^'^'*' 
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and A*^'^^^ are present, while for mixed neutron-star- 
black-hole systems only one term, corresponding to the 
neutron star, is present; see below]. Here, we consider 
a binary system of (gravitational) masses niA and m^, 
and denote 



AI = rriA + ms 



rriA ruB 
{niA + msY 



(1.2) 



[A labelling of the two bodies by the letters A and B will 
be used in this Introduction for writing general formulas. 
We shall later use the alternative labelling A = 1, B = 2 
when explicitly dealing with the metric generated by the 
two bodies.] The binary black-hole (or point mass) po- 
tential A^^^(r) is known up to the third post-Newtonian 
(3PN) level [1^, namely 

A^p^ir) = l-2u + 2i^u^ + aii^u'^ , (1.3) 

where 04 = 94/3 - (41/32) vr^ ~ 18.68790269, and 



GM 



(1.4) 



shall always put the multipolar index £ within paren- 
theses to avoid ambiguity with our later use of the la- 
belling ^, _B = 1, 2 for the two bodies.] The correspond- 
ing leading-order radial potential of the gravito-magnetic 
type is proportional to it^^+'^ (instead of u^^"*"^), and 
to j'^^ R^~^^, where denotes a dimensionless "mag- 
netic tidal Love number" . It was found 0, 0] that both 
types of Love numbers have a strong dependence upon 
the compactness Ca = G niA/ {c^ Ra) of the tidally de- 
formed body, and that both k'^^ and j^-* contain a factor 
l — 2CA^ so that they would formally vanish in the limit 
where body A becomes as compact as a black hole (i.e. 
Ca Cbh = \)- This is consistent with the decompo- 
sition (11.11) . where the binary black- hole radial potential 
A^^^{r) is the only remaining contribution when one 
formally takes the limit where both Ca and Cb tend to 
the black- hole value Cbh = 1/2. Finally, the supplemen- 
tary factors A^^iieetric('^) ^nd li^Lg„etic(^) dcuote the 
distance-dependent amplification factors of the leading- 
order tidal interaction by higher-order PN effects. They 
have the general form 



with c being the speed of light in vacuum and G the 
Newtonian constant of grav itation. 

It was recently found |22l . [23j that an excellent descrip- 
tion of the dynamics of binary black-hole systems is ob- 
tained by augmenting the 3PN expansion Eq. (jl.3p with 
additional 4PN and 5PN terms, and by Pade resumming 
the corresponding 5PN Taylor expansion. 

The tidal contributions A^j^^{r) can be decomposed 
according to multipolar order and type, as 



A 



tidal 



ir)=Y.[^'2.\L.ir)A^^ 



+ A 



J^W LO 

l>2 
(£)LO 



A magnetic V ) A magnetic 



lectricV / 



}• (1-5) 



Here, the label "electric" refers to the gravito-electric 
tidal polarization induced in body A by the tidal field 
generated by its companion, while the label "magnetic" 
refers to a corresponding gravito-magnetic tidal polar- 
ization. On the other hand, the label LO refers to the 
leading-order approximation (in powers of u) of each 
(electric or magnetic) multipolar radial potential. For 
instance, the gravito-electric contribution at multipolar 
order I is equal to @ 



electric - " 



where 



■V A — ^ ru A 



niA \G{mA + tub) 



21+1 



(1.6) 



(1.7) 



Here, Ra denotes the radius of body A, and fc^' de- 
notes a dimensionless "tidal Love number". [Note that 

k'^P was denoted kf in our previous work. Here we 



^Velectric(^) — 1 + Q^lelictric" + electric'^"' + . . . , (1 



Aii) 



Aie) 



(r) 

A magnetic V / 



1 



A{1) 
^1 magnetic^ 



(1.9) 



where u is defined by Eq. (|1.4I) above. 

The main aim of the present investigation will be to 

compute the electric- type amplification factors electric' 
for ^ = 2 (quadrupolar tide) and £ = 3 (octupolar tide), 

at the second order in u, i.e. to compute both af^fl^^^.^^ 



and 



*2elcctr 



We shall also compute the magnetic-type 

amplification factor magnetic for ^ = 2, at the first 
order in u. 

The analytical value of the first-order electric amplifi- 
cation coefficient ctfl,fl^^- was computed some time ago 
for ^ = 2 (see Ref. [29] in fl]) and was reported in Eq. (38) 
of 3, namely 



Aie=2) _ 5 

"l electric " 2 



(1.10) 



where Xa = mA/{mA + uib) is the mass fraction of 
body A. The analytical result (|1.10l) has been recently 
confirmed @. On the other hand, several comparisons of 
the analytical description of tidal effects with the results 
of numerical simulations have indicated that the ampli- 



fication factor A 

A{i=2) 



e=2) 

A electric 



[r) is larger that its IPN value 
1 -|- a^li^t^ic ^^'^ have suggested that the higher-order 
take large, positive values. More 



coefficients a 



A{i) 

2 electric ' " 



precisely, the analysis of Ref. [5] suggested (when taking 
into account the value (jl.lOp for ai) a value of order 

'^2^ofe^tric ^ + 40 (for the equal-mass case) from a com- 
parison with the numerical results of Refs. [l^, [2l| on 
quasi-equilibrium adiabatic sequences of binary neutron 



3 



stars. Recently, a comparison with dynamical simula- 
tions of inspiralling binary neutron stars confirmed the 
need for such a large value of a2^eiectric I, [Note 
that, while the comparison to the highest resolution nu- 
merical data suggests the need of even larger values of 



A(i=2) 

V 

2 electric 



of order + 100, the comparison to approximate 
resolution-extrapolated data only call for 012 values of or- 
der + A0. See Fig. 13 in [10|].] 



II. EFFECTIVE ACTION APPROACH TO 
TIDAL EFFECTS 

It was shown long ago 24], by using the technique 
of matched asymptotic expansions, that the motion and 
radiation of N (non-spinning) compact objects can be 
described, up to the fifth post-Newtonian (5PN) approx- 
imation, by an effective action of the type 



5*0 



IGttG 



s/gR{g) + Spoint mass ; 

(2.1) 



where R{g) represents the scalar curvature associated 
with the metric g^v, with determinant —g, and where 



point mass 



E 



niA (? dTA 



(2.2) 



is the leading order skeletonized description of the com- 
pact objects, as point masses. Here dxA denotes the 
proper time along the worldline UAiTA) of A, namely 
dTA = c^^{-gf^i,{yA)dy'^dy'^)^/'^. To give meaning to 
the notion of point mass sources in General Relativity 
one needs to use a covariant regularization method. The 
most convenient one is dimensional regularization, i.e. 
analytic continuation in the value of the spacetime di- 
mension D = 4 -\- s, with s € C being continued to zero 
only at the end of the calculation. The consistency and 
efficiency of this method has been shown in the calcula- 
tions of the motion [2^ [2^ and radiation of binary 
black holes at the third post-Newtonian (3PN) approxi- 
mation. 

It was also pointed out in [23] that finite-size effects 
(linked to tidal effects, and the fact that neutron stars 
have, contrary to black holes, non-zero Love numbers 
fc^-*) enter at the 5PN level. In effective field theory, 
finite-size effects are treated by augmenting the point- 
mass action (j2.2p by nonminimal worldline c oupling s in- 
volving higher-order derivatives of the field [28l - l3(j |. In 
a gravitational context this means considering worldline 
couplings involving the 4- velocity = dy'^ /dTA (satisfy- 



ing g^u "a '^a 



'C ) together with the Riemann tensor 



RaPiiv and its covariant derivatives. To classify the pos- 
sible worldline scalars that can be constructed one can 
appeal to the relativistic theory of tidal expansions [sil - 
l33j . In the notation of Refs. [33, [s^ the tidal expansion 
of the "external metric" felt by body A can be entirely 
expressed in terms of two types of external tidal gradi- 
ents evaluated along the central worldline of this body: 



the gravito-electric Gf (r^i) 



Ga,...ai{'rA) and gravito- 
magnetic H^{ta) = ^^{ta) symmetric trace-free 
(spatial) tensors, together with their time-derivatives. 
[The spatial indices a.i = 1,2,3 refer to a local frame 
X\ = CTA, X\ attached to body A] This implies that 
the most general nonminimal worldline action has the 
form 



'5'nonminimal — 



Y.Y.{\1^^'^ IdrAiGtiTA)? 

a'i^J dTA{Ht[TA)f 



A £>2 

1 I I I 

2 l+l J\ ~^ 

1 1 1 ,/W 

2 i\ c 



^T^7l^'A' I dTA{Gt{TA)r 



I e 1 I 



1(1) 



2 £+1 £1 



^A 



dTA{Ht{TA)) 



(2.3) 



where Gf (r^) = dG^/dTA, and where the ellipsis refer 
either to higher proper-time derivatives of Gf^ and H^, or 
to higher-than-quadratic invariant monomials made from 
Gf , and their proper-time derivatives. For instance, 
the leading-order non-quadratic term would be 



dTA G'^ij G^ G^ 



(2.4) 



Note that the allowed monomials in Gl, Hl and their 
time derivatives are restricted by symmetry constraints. 
When considering a non-spinning neutron star (which is 
symmetric under time and space reflections) one should 
only allow monomials invariant under time and space re- 
versals. For instance Gab Gab and Gab Hab are not al- 
lowed. 

The electric-type tidal moments G^ are normalized 
in a Newtonian way, i.e. such that, in lowest FN or- 
der, they reduce to the usual Newtonian tidal gradients: 
Gf = [dL C/(X")]x»=o + (^), where U{X) is the New- 
tonian potential and = di^ . . . di^ represents multi- 
ple ordinary space derivatives. The magnetic-type ones 
Hf^ are defined (in lowest FN order) as repeated gradi- 
ents of the gravitomagnetic field c^goa- With these nor- 



malizations the coefficients /^]^ and 



(T^i^ in the nonmin- 



imal action (|2.3p both have dimensions [length] ^^"'"^/G. 
They are related to the dimensionless Love numbers k 
and j^-*, and to the radius of body A, via [^ 



{2£- 



1)!! 



R 



2i+l 
A ■ 



Ga 



1 



4(^4-2) (2£- 1)!! 



]a ^/ 



(2.5) 



(2.6) 



Note that the coefficients associated with the first time 
derivatives of Gf and have dimensions G^'^^ ~ 



4 



[length] 



2f+3 



Ga 



The nonminimal action (12.31) has 



a double ordering in powers of Ra and in powers of 1/c^ . 
The lowest-order terms in the Ra expansion are propor- 
tional to R\ and correspond to the electric and magnetic 
quadrupolar tides, as measured by and i?^, respec- 
tively. 

We have written the most general nonminimal action 
(|2.3p in terms of the irreducible symmetric trace- free spa- 
tial tensors (with respect to the local space associated 
with the worldline ?/^(t^)) describing the tidal expan- 
sion of the "external metric" felt by body A, as defined 
in [s^]- These tidal tensors played a useful role in sim- 
plifying the (IPN- accurate) relativistic theory of tidal ef- 
fects. In our present investigation, it will be convenient 
to express them in terms of the Riemann tensor and its 
covariant derivatives. Eq. (3.40) in [s^l shows that (in 
the case where one can neglect corrections proportional 
to the covariant acceleration of the worldline) the first 
two electric spatial tidal tensors. Gab and Gabc, are sim- 
ply equal (modulo a sign) to the non-vanishing spatial 
components (in the local frame) of the following space- 
time tensors (evaluated along the considered worldline) 



G 



af3 



R. 



GajB-t = - SymQ,^^(V^ Ri3fj.f,^) 



(2.7) 



(2., 



Here the notation Ga/B for (minus) the electric part 
of the curvature tensor should not be confused with 
the Einstein tensor, Sym^^^^ denotes a symmetrization 
(with weight one) over the indices a/37, while V;^ = 
P(u)'^ denotes the projection of the spacetime gradi- 

-2 



ent orthogonally to u'^ {P{u)'l = 
[Note that in the Newtonian limit u^^ ~ cSq so that 
the Newtonian limit of Gap is — c^i?ao/3o, where the 
factor cancels the 0(l/c^) order of the curvature 
tensor.] By contrast, the presence of the extra term 
- 3 E*^^ E*^ on the right-hand side of Eq. (3.40) in 
[33 | shows that the £ = 4 electric spatial tidal tensor 
Gabcd = di^abc E'^) would differ from the symmetrized spa- 
tial projection of (VqV^ R-yfj.Su) u^^ by a term propor- 
tional to Gi^a-y Gps) . [Here, the angular brackets denote 
a (spatial) symmetric trace-free projection.] In addition, 
the electric time derivatives, such as Gab can be replaced 
by corresponding spacetime tensors such as Gap ■ 

Similarly to Eqs. (UTT]), (H^, one finds that the £ = 2 and 
^ = 3 magnetic tidal tensors (as defined in Refs. [3^l33|) 
are equal to the nonvanishing local-frame spatial compo- 
nents of the spacetime tensors 



Hap^+2cR, 



(2.9) 



Hap^ = + 2 c Sym„^^(Vi R}^^,) u'^ u'' , (2.10) 

where R*^^ai3 = 5 ^^^:^po■ -R'"^ ap is the dual of the curva- 
ture tensor, e^jypo- denoting here the Levi-Civita tensor 



(with 60123 = +^/g)■ Note the factor +2 entering the 
link between the magnetic tidal tensors Hap , ■ ■ ■ (nor- 
malized as in Refs. [32, 113]) and the dual of the curva- 
ture tensor, which contrasts with the factor —1 entering 
the corresponding electric tidal-tensor links, Eqs. (|2.7p . 
([2:81) . [The definition of B^^ in the text below Eq. (5) of 
should have included such a factor 2 in its right-hand 
side. On the other hand, the corresponding magnetic- 
quadrupole tidal action, Eq. (13) there, was computed 
with Hab and was correctly normalized.] Let us also note 
that the expressions (|2.7p - (|2.10p assume that the Ricci 
tensor vanishes (e.g. to ensure the tracelessness of Gap)- 
One could have, alternatively, defined Gap etc. by using 
the Weyl tensor Gafipu instead of RafiPv However, as 
discussed in pol j. the terms in an effective action which 
are proportional to the (unperturbed) equations of mo- 
tion (such as Ricci terms) can be eliminated (modulo 
contact terms) by suitable field redefinitions. 

Finally, the covariant form of the effective action de- 
scribing tidal interactions reads 

Stot -\- 'S'point mass ^" '5'i;ionminimal (^'i-i-) 

where Sq and S'pointmass are given by Eqs. (|2.ip . (12.21) . 
and where the covariant form of the nonminimal world- 
line couplings starts as 



1 



>5'nonminimal — ^ ^ { T /^A ^ / '^^A 
A ^ 

A fjaP 
A 



^A /-iOtP 



-^2^A j dTAH^pH 
j dTAGtp.G 
4''! dTA{u^^V,Gtp){u\VMf) 



aP^ 
A 



J_ ,/(2) 

4c2 



(2.12) 



where G'^ = 9°''^ g^'^ G^^, etc. [evaluated along the A 
worldline] . 

In principle, one can then derive the infiuence of tidal 
interaction on the motion and radiation of binary sys- 
tems by solving the equations of motion following from 
the action (|2.1ip . (|2.12p . More precisely, this action 
implies both a dynamics for the worldlines where the 
geodesic equation is modified by tidal forces (coming 
from (55nonminimai/^ 2/^^('''a)), and modified Einstein equa- 
tions for the gravitational field of the type 



Ru 



1 



R gi-iu 



SttG 



^^point 1 



T, 



nonminimal! 



where the new tidal sources 

^'nonminimal 



nonminimal V / 



J 1 

(2.13) 



(2c/7g)(SS'"°"'°'"™'^755^^(a;) are, essentially, sums 
of derivatives of worldline Dirac-distributions: 



T, 



nonminimal 



A I 
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The task of solving the coupled dynamics of the world- 
lines and of the gravitational field, both being modified 
by tidal effects, at the 2PN level, i.e. at the next-to- 
next-to-leading order in tidal effects, and then of com- 
puting the looked for higher-order terms in the amplifi- 
cation factors ()1.8p . ()1.9p is quite non-trivial. Happily, 
one can drastically simplify the needed work by using a 
general property of reduced actions. Indeed, we are in- 
terested here in knowing the influence of tidal effects on 
the reduced dynamics of a compact binary, that is the 
dynamics of the two worldlines ^^(t), ^^(t), obtained 
after having "integrated out" the gravitational field (i.e. 
after having explicitly solved g^i, (x) as a functional of the 
two worldlines). When considering, as we do here, the 
conservative dynamics of the system (without radiation 
reaction), it can be obtained from a reduced action, which 
is traditionally called the "Fokker action". See [2^ and 
references therein for a detailed discussion (using a dia- 
grammatic approach) of Fokker actions (at the 2PN level, 
and with the inclusion of scalar couplings in addition to 
the pure Einsteinian tensor couplings). If we denote the 
fields mediating the interaction between the worldlines 
y = {yA,yB} as ip (in our case cp = g^v), the reduced 
worldline action S'red[y] (a functional of the worldlines y) 
that corresponds to the complete action S[ip,y] describ- 
ing the coupled dynamics of y and Lp is formally defined 
as: 

5,.ed[y] = 5[^soi[2/],2/], (2.14) 

where (psoi [y\ is the functional of y obtained by solving 
the (/3-field equation, 

5S[vM/5v^^. (2.15) 

considered as an equation for tp, with given source- 
worldlines. [This must be done with time-symmetric 
boundary conditions, and, in the case of g^^, the addition 
of a suitable gauge-fixing term; see [l^ for details.] 

Having recalled the concept of reduced (or Fokker) ac- 
tion, let us now consider the case where the complete 
action is of the form 

5[vp,y] = + (2.16) 

where e denotes a "small parameter" . In our case, e 
can be either a formal parameter associated with all 
the nonminimal tidal terms in S'nonminimai, Eq. (|2.12p . 
or, more concretely, any of the tidal parameters entering 
(|2.12p : /^^~^', A'b~^\ 6tc. As said above, when turning 
on e, the equations of motion, and therefore the solu- 
tions of both (/J and y get perturbed by terms of order e : 
ip = (^(0) -(-g (^(1) -(- . . y = y(o) _(_g -(-... but a simplifi- 
cation occurs when considering the reduced action (|2.14p . 
Indeed, it is true that the field equation (|2.15|) for gets 
modified into 

^_ 5S[pM 5S(^)[p.,y] ^ ^ 5S''^)[pM ^ ^^.ij) 
5 ip 5p} 5 ip ^ 



so that its solution Lp^oi [y\ gets perturbed: 

Vsox[y]^p'^^[y]+ep^^M+0{e'). (2.18) 

However, when inserting the perturbed solution (|2.18p 
into the complete (perturbed) action (|2.16p one finds 

5,edM = S[p^M + ^v's^M + 0{e^),y] 

= S [yU] + e pi'2 [y] f| ['Piol [yly] + 0{e') 
= S[p^^,[yly] 

+^^iol[y]^[^lol [yU] + o{e') 

= S[p^^,[yly] + Oie'), (2.19) 

because, by definition, (^^^J is a solution of (5 5*^"-' /S p ~ 0. 
The final result is therefore 

5„d[y] = 5Wb(::|M,y] + e5W[^(°|M,2/] + 0(e2) 
= Si'j^ [y] + eS('^[p^^2 [yly] + 0(6^) . (2.20) 

In words: the order 0{e) perturbation 

eSll^,[y]^S,M~Si%] 

of the reduced action is correctly obtained, modulo terms 
of order O(e^), by replacing in the 0(e) perturbation 

e^W [p,y] 

of the complete (unreduced) action the field p by its un- 
perturbed solution yjgl^j [y] . 

In our case, the ordering parameter e is either the col- 

1 (2) (2) (3) (3) (2) _2 '(2) _2 

section , fi'g\ fiy , fj,)^\ ay c c ^, . . ., 

or the corresponding sequence of powers of Ra and 
Rb ■ R'a, R%, R\, ■ • . The terms quadratic in e would 
therefore involve at least ten powers of the radii (and 
would mix with higher-than-quadratic worldline contri- 
butions akin to (|2.4p '). Neglecting such terms, we con- 
clude that the higher-PN corrections to the tidal effects 
are correctly obtained by replacing in Eq. (|2.12p . con- 
sidered cLS a functional of {x) and y^(r^), the metric 
gtiu{x) by the point-mass metric obtained by solving Ein- 
stein's equations with point-mass sources. [This was the 
method used by one of us (T.D.) to compute the IPN 
coefficient (|1.10p from the calculation by Damour, Soffel 
and Xu of the IPN-accurate value of Gab [s^. Is^.] 



III. THE 2PN POINT-MASS METRIC AND ITS 
REGULARIZATION 

The result of the last Section allows one to compute 
the tidal corrections to the reduced action for two tidally 
interacting bodies A, B with the same accuracy at which 
one knows the metric generated by two (structureless) 
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point masses inA,yA:''^B,yB- The metric generated by 
two point masses has been the topic of many works over 
many years. It has been known (in various forms and 
gauges) at the 2PN approximation for a long time [36'- 
|38|. Here, we shaU use the convenient, exphcit harmonic- 
gauge form of Rcf. [39], with respect to the (harmonic) 
coordinates a;^ = (x*^ = ct. a;*), i.e. the metric 

ds'^ = goo{dx")'^ + 2 goi dx^dx' + g^-, dx'dx' , (3.1) 

where, at 2PN, the metric components are written as 

500 = -l + 2e^V -2e^V^ 

' 0(8), 



9oi 



Ae^Vi-8e'^R, + 0{7), 



g.^j = S,j {l + 2e^V + 2e^V'^) +Ae^ W^j + 0(6) . 

(3.2) 

Here, as below, we sometimes use the alternative notation 
e = 1/c for the small PN parameter. We used also the 
shorthand notation 0{n) = 0(e") = 0(c-"). 

The various 2PN brick potentials V, Vi , Wij , Ri and X 
are the (time-symmetric) solutions of 



4 TT Gcr, , 



4TrG{Va^-V,a)- 2 dkV O^Vk - - dtV O^V , 



ax ^ ~A^TGVa^^ + 2V^^t^^V + Vd^V +-{dtVf 

(3.3) 



where dt denotes a time derivative (while we remind that 
di, for instance, denotes a spatial one), and where the 
compact-supported source terms are [40| 

2^00 I r^ii rpQi 

<J= 2 , a^ = , a^J=r^ (3.4) 

with T'^" being the stress-energy tensor of two point 
masses: 

T^-- = ^il{t) <(i) <(t) 5{x - yi{t)) + 1^2, (3.5) 



where 



/ii(t) = mi \g-^'\g^uv^^v'i/^)-'/^ 



(3.6) 



Here, = = {c,v\) and the index 1 on the bracket 
in Eq. (13.61) refers to a regularized limit where the field 
point tends towards the (point-mass) source point y\. 
Note that, in this section, we shall generally label the 
two particles as {mi,y\), {1112, y 2), instead of (m^,?/^), 
{mB,yB) as above. The notation 1 o 2 means adding 
the terms obtained by exchanging the particle labels 1 
and 2. 

The explicit forms of the (time-symmetric) 2PN- 
accurate brick potentials V, Vi, etc. are [11] 



V = ^ + ^(J-!^+^-^ + Gm2(-'' 
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(3.9) 



6 _ ^ / ("12^1 ) 1 , M 2(n2^;i) 3(n2t^1 

i?, - G 777 17^277^2 I ^ (^^J + ^ J " + | 

+ -1 {^^^fpl + (2(7Zi2.l) - + 2(7.2.1) - ^) } 

+ A + G^mi7772 + ^) I - -^^^ + 0(2) + 1^2, (3.10) 



^ = ^ ((-1-0^ - + G^mi7n2.? + + 
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+ 0(2) + 1 o 2 . 



Here ri = x - yi, ri = \ri\, rii = ri/ri, r2 = x - y2, 
etc., yi2 = yi - 2/2, 7-12 = |yi2|, "12 = yi2/ri2, vi2 = 
i>i — V2, (ni2 vi) = ni2 ■ Vi, while the quantity S denotes 
the perimeter of the triangle defined by x, yi and y2, viz 

S = ri+r2+ri2. (3.12) 

In all the PN expressions, the spacetime points x^,yi, yi^ 
(and the velocities v^) are taken at the same instant 
i.e. x'^ = yi = y2 = ct. 

Let us now discuss in more detail the crucial operation 
(already implicit in Section |TT] above) of regularization 
of all the needed field quantities, such as g^i,{x), g{x), 
Rfiai^isix), . . ., when they are to be evaluated on a world- 
line: — > 2/^. As mentionned at the beginning of Sec- 
tion all the quantities [G^i/(a;)]i, . . . , [Rfiai^i3{x)]i are 
defined by dimensional continuation. It has been shown 
long ago 0,|4l| that, at 2PN, dimensional regularization 
is equivalent to the Riesz' analytic regularization, and is 
a technical shortcut for computing the physical answer 
obtained by the matching of asymptotic expansions. In 
addition, because of the restricted type of singular terms 
that appear at 2PN (see Eqs. (25), (30) and (33) in 
the analytic-continuation regularization turns out to be 
equivalent to Hadamard regularization (used, at 2PN, in 
[H, [3^ 113)) see below. Here, it will be technically con- 
venient to use Hadamard regularization (which is defined 
in D = 4) because the explicit form p.7p - (l3.1ip of the 
2PN metric that we shall use applies only in the physical 
dimension D = 4 and has lost the information about its 
dimensionally continued kin in _D = 4 + e. 

Let us summarize here the (Hadamard-type) defini- 
tion of the regular part of any field quantity ip{x) (which 
might be a brick potential, V{x), Vi{x), . . a component 
of the metric g^^{x), or a specific contribution to a tidal 
moment. Gap,---)- We consider the behavior of ip{x) 
near particle 1, i.e. when ri = |x — yi | — >■ 0. To ease 
the notation, we shall provisionally put the origin of the 
(harmonic) coordinate system at yi (at some instant t), 
i.e. we shall assume that yi — 0, so that ri = \x\ = r 
and rii = ri/ri = x/r = n. We consider the expan- 
sion of '-p{x) in (positive and negative) integer powers k 
of ri — 7', and in spherical harmonics of the direction 
ni = n, say (for /c e Z, ^ € N, TV € N) 

k>-N l.>0 

where = n"-i---°-t denotes the symmetric trace- free pro- 
jection of the tensor = n"-^ . . . n"*. [The angular func- 

+1 

tion /I; is equivalent to a sum of ^ c„i Yi„i ■] We 

m=-£ 

(uniquely) decompose the field ^p{x) in a regular pait (R) 
and a singular one (S), 

ifix) = R[cp{x)]+S[ip{x)], (3.14) 



(3.11) 

I 

by defining {n £ N) 

i?Ka.)]EE^^/+2"n^/^+2", (3.15) 

l>0 n>0 

SMx)]^ ^ r'^n^fl (3.16) 

Note that R [(p{x)] can be rewritten as a sum of infinitely 
differentiable terms of the type x^{x'^)^. By contrast 
S[(p{x)] is such that it (if N, in Eq. (IXTO)) . is strictly 
positive), or, one of its (repeated) spatial derivatives, 
tends towards infinity as r — >■ 0. Note also that the 
R + S decomposition commutes with linear combinations 
(with constant coefficients), as well as with spatial deriva- 
tives, in the sense that R [a (p{x) + bil;{x)] = aR [(p{x)] + 
bR[i;{x)], S[aLp{x) + bi}{x)] ^ a S [ip{x)] + b S [ij{x)], 
R[di(p{x)] — diR[tp{x)] and Sldif^x)] = diS[(p{x)]. 
By contrast, the R + S decomposition (as defined above, 
in the Hadamard way) does not commute with nonlin- 
ear operations (e.g. R[(pip] 7^ R[ip] R[iIj]), nor even 
with multiplication by a smooth (C°°) function f{x) (e.g. 
R [f(p] / i? [(/?]). This is a well-known inconsistency of 
the Hadamard regularization, which created many ambi- 
guities when it was used at the 3PN level [i^ |3] . One 
might worry that our present calculation (which aims at 
regularizing nonlinear quantities quadratic in R^ai^p ^ 
9^9 + dgdg) might display such ambiguities already 
at the 2PN level. Actually, this turns out not to be the 
case because of the special structure of the 2PN met- 
ric which is at work in the Riesz-analytic-continuation 
derivation of the 2PN dynamics in [2J| . Let us first 
recall why the Riesz-analytic-continuation method, or, 
equivalently (when considering the regularization of the 
metric and its derivatives), the dimensional-continuation 
method, is consistent under nonlinear operations. The 
dimensional-continuation analog of (I3.14p - p.l6p consists 
of distinguishing, within ^{x), the terms that (in dimen- 
sion 4 + e) contain powers of r of the type r'^~"^, with 
n = 1, 2, 3, . . . (which define the ^-singular part of ^p{x)), 
and the terms that are (formally) C°° in 4 -I- £ dimen- 
sions (which define the e-regular part of ip{x)). It is then 
easily seen in dimensional continuation (simply by con- 
sidering the continuation to large, negative values of the 
real part of e) that the e-singular terms give vanishing 
contributions when evaluated at r — 0, and that they 
do so consistently in nonlinear terms such as dip dtp. Let 
us now indicate why the special structure of the 2PN 
metric ensures that the decomposition into e-singular 
parts and e-regular parts of the various brick potentials 
V{x), Vi{x), . - . coincides with their above-defined decom- 
position into Hadamard-singular [S \V{x)], S \Vi{x)], . . .) 
and Hadamard-regular parts {R\V{x)\,R\Vi{x)\, - . -) in 
the four-dimensional case. This is trivially seen to be 
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the case for most of the 2PN contributions to the brick 
potentials (because one easily sees how those contribu- 
tions smoothly evolve when analytically continuing the 
dimension). However, the most nonlinear contributions 
to the 2PN metric, namely the terms, say X^^^^\ in X 
that are generated by the cubically nonlinear terms con- 
tained in the last source term, W^^^ dijV , on the right 
hand-side of the last Eq. p.3p (where W-^"^^ is the part 
of Wij generated by —diVdjV) are more dehcate to dis- 
cuss. Actually, among the contribution X^^^^\ only the 
terms proportional either to m\m2 or to toito|, i.e. the 
terms whose cubically nonlinear source ~ d^Vl^~^dVdV 
involve two V potentials generated by one worldline and 
one V potential generated by the other worldline, such as 
(X m\m2, pose a somewhat delicate problem. 
More precisely, it is easily seen that the only dangerous 
part in 

X^ViViV^)^ considered near the first worldline, is 
of the form f [x) / r^^^"^^^ in dimension 4 -I- e, where fix) 
denotes a smooth function. [Here, we add back the parti- 
cle label indicating whether the expansions p.lSp . p.l6p 
refer to the first {A = 1), or the second {A = 2) particle. 
The appropriate label should be added both on r and n 
in Eqs. ^Jb^i, ([XTC)) : r'' r\h\] The problem is 
that the power of 1/ri in this e-singular term becomes 
an even integer when e — ;> 0. When inserting the Taylor 
expansion of f{x), say f{x) ^ I] /^+^", some 

of the terms in the e-singular contribution f [x] / r^^'^'^^^ 
might be of the form r^+2"'-2e ^^--^^y^ n' = n - 1 > 0, 
and might then contribute to the Hadamard-regular part 
of in the hmit e 0.^ This would mean that 

the Hadamard-regular part of x'--^^^^^'^) would not co- 
incide with its e-regular part. We already know from 
Refs. [H, d^), which used Hadamard regularization to 
derive the 2PN-accurate dynamics and found the same 
result (modulo gauge effects) as the analytic-continuation 
derivation of Ref. |2J], that this is not the case for the 
regularized values of X^^^^^^^^ and of its first derivatives 
on the first worldline. [Indeed, these quantities enter the 
computation of the equations of motion.] On the other 
hand, the computations that we shall do here involve 
higher spatial derivatives of X, and it is important to 
check that we can safely use Hadamard regularization to 
evaluate them. This can be proven by using the tech- 
niques explained in [24] . based on iteratively consider- 
ing the singular terms in W^-J^"^ and X^^^^^ generated 
by the singular local behaviour (near the first world- 
line) of their respective source terms. One finds then 
that the smooth function f(x) entering the dangerous 

terms f Ix) / r^^^"^^^ in X^^^^^^^^ is of the special form 
f{x) ^ J2^i^Lrinf in dimension 4 + e, with £ > 1, 
where Gl = 9lV2 denotes the £-th tidal gradient (consid- 
ered near the first worldline) of the V potential generated 
by the second worldline. When working (as we do) at the 
2PN accuracy, we can take V at Newtonian order, and 
the gradients Gl — [dL{Gm2/r'^~^^'')]i are then trace- 
less: Gl = Ga^a2 -ai = G {ata^-ai) ■ As a couscquence. 



it is immediately seen that, in the limit e — )■ 0, the po- 
tentially dangerous term /(a;)/rp^^^^ in X^^'^-^'^^^) does 
not give any contribution to the Hadamard-regular part 
of X. This means that we can compute the £-regularized 
reduced tidal action (|2.12p by replacing, from the start, 
the brick potentials V{x)^ Vi{x), . . by their Hadamard- 
regularized counterparts, R \V{x)], R \Vi{x)], . . . 

Summarizing: The A-worldline part of the tidal ac- 
tion (|2.12p can be obtained by computing all its elements 
[dTA = c-^{~gf,,{yA) dy'Xdy'Xy^^ ,G^p, . . .) within the 
A- regular metric g'^^J''^'^^ {x) obtained by replacing each 
2PN brick potential V{x), Vi{x), ... by its A-Hadamard- 
regular part Ra[V{x)], RA[Vi{x)], . . . 

As a check on our results we have also computed the 
reduced action by effecting the Hadamard regularization 
only at the end, i.e. by computing, e.g., RaIGop G"^], in- 
stead of {Gafi G°'^)[RA{V),RA{Vi), . . .]. Note that this 
second calculation involves a more nonlinear mixing of 
Hadamard-regular and Hadamard-singular parts than 
the first one. However, as the Hadamard-delicate terms, 
at 2PN, in G^/? G"^ are of the type d^jVd^jX'-^^^^^^\ 
one can check (using the special structure of the danger- 
ous terms in given above) that effecting the 
Hadamard regularization at the end will agree with the 
first method where the Hadamard regularization is done 
on the brick potentials. This result hods true also for the 
higher multipolar tidal-action terms such as Gap^y G"^^ . 
[Note in passing that this might no longer be the case 
if one were considering the regularization of a term cu- 
bic in the tidal potentials, such as Gap G^'^G'".] In the 
present case, our final results were indeed found to be the 
same, thereby checking the many complicated algebraic 
operations necessary to derive them. 

Let us now give some indications on the computa- 
tion of the regular parts of the various brick potentials 
V{x), Vi{x), . . . The situation is very simple for the "lin- 
ear potentials" V and Vi, which satisfy linear equations 
with delta-function sources (see Eqs. p.3p ). Near, say, 
the particle A = 1, the A-regular parts of V and Vi 
are the terms in Eqs. p.7p . p.8p which are generated 
by the source terms oc 6{x — 1/2) of the second particle. 
It is indeed easily seen (from the definition p.lSp ) that 
the 1-regular part of all the terms explicitly written in 
Eq. p.7p vanishes, while all the non-explicitly written 
terms obtained by the 1^2 exchange are regular near 
the particle 1. The same is true for Vi, Eq. p.Sp . A sim- 
ple rule for obtaining these results is to note that, from 
the definition p.l6p . any term of the form 

rf+V(a;), fceZ, (3.17) 

where f{x) is a smooth function of x'^ (near x — yi 
at fixed instant t), and where the power of ri is odd, is 
purely singular. 

The situation is more complicated for the higher-order 
potentials Wij and Ri, whose sources contain both com- 
pact terms oc S{x — yA), and quadratically nonlinear non- 
compact ones cx dV dV, and still more complicated for 
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the X potential whose source even depends on the pre- 
vious Wij potential. 

The potential Wij can be decomposed in powers of the 
masses. It contains terms proportional to mi, TO2, itt-i, 
and mi m2 . It is easily seen that while the terms propor- 
tional to TOi and mf are 1-singular, the terms propor- 
tional to TO2 and m| are 1-regular. It is more delicate to 
decompose the mixed terms oc mi m2 into 1-regular (i?i ) 
and 1-singular (Si) parts. More precisely the mi m2 part 
of Wij has the form 



it is 1-regular, but might be singular near particle 2.] Be- 
cause S — ri +r2 + ri2 is a function of "mixed character" 
(partly regular and partly singular), it is not immedi- 
ate to decompose the functions p. 221) into 1-regular and 
1-singular parts. [This mixed character of S is deeply 
linked with the fact that it enters the 2PN metric be- 
cause of the basic fact that a solution of Ag = is 
g = \tl S.] a simple (though somewhat brute force) way 
of extracting the regular parts of the functions p.22p 
consists of decomposing S into 



■u(o) 



ij(0) 



S = So + ri = Sq { 1 



So 



(3.23) 



where 



with 
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\o ri2S 
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ri2S 



P{ni2) 



J_ (nfni^ + 2nl'n?2 - n{2n{2) , (3.19) 



' «J (0) 



r^S [S^ 

^^12^12 1^^ + 
1 



ri2S 



P{ni2r 



+ ^ [^2^^"! - 2^ «i2 - <2"i2) , (3.20) 



and where P(ni2)* 



5^] 



'12 "-12 



denotes the projec- 



tor orthogonal to the unit vector ni2. [The decomposi- 
tion (I3.18P simply corresponds to the decomposition of 
Eq. p. 101) into an explicitly written term and its 1 O 2 
counterpart.] Here we see that there appear (modulo 



x-independent factors, such as r 
terms of the type 
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12 I 
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52 
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^2 



52 



P(ni2)^-',...) 



(3.21) 



where we recall that S = ri + r2 + ri2- Near particle 1, 
rij is a smooth function, while n\ = r\/ri is the ratio of 
a smooth function {r\ = a;* — j/i) by ri. In other words, 
the five terms listed in Eq. p.2ip are of three different 
types: 
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52 



and 
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S2 ' 



(3.22) 



5*0 = r2 + ri2 



(3.24) 



(note that So is a smooth function near particle 1), and 
then expanding S"^" in powers of ri/ So- Namely 



-=-(i--+4-4 

S So \ So Sq Sq 
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(3.26) 



and more generally 
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(3.27) 



Using these expansions, together with the rule that terms 
of the form p.l7p are purely singular, it is easy to derive 
the following results for the 1-regular parts of functions of 
the type p. 221) . and, more generally, of the types f{x)/S, 
f{x)/S\ f{x)/{r, S) and f{x)/{n S'): 
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(3.28) 



where f{x) denotes a generic smooth function near par- 
ticle 1. [As we always consider the neighbourhood of 
particle 1 we do not add an index to f{x) to recall that 
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Here we use a lower R subscript ((p(x))ii to denote 
the 1-regular part of a function ip{x) (above denoted as 
R[lp{x)]). [We omit to decorate R with a label 1, but 
one should remember we are always talking about the 
1-regular part of f{x).] 

Note that, as indicated, all the terms above have the 
simple property that the regular-projection operator R 
commutes with the multiplication by a smooth function, 
e.g. R \f{x) 5'-!] = f{x) R Beware that this prop- 

erty is only true for the special singular terms considered 
here. We shall later see that more complicated singu- 
lar terms (entering the X potential) do not satisfy this 
simple commutativity property. 

Note that the number of terms one needs to retain in 
the above expansions depends on the quantity one wants 
to evaluate on the first worldline. For instance, when 
evaluating G^^, which involves the curvature tensor, and 
therefore two spatial derivatives of the metric (and, in 
particular, of i?[Wij]) we need to include enough terms 
to ensure that i?[Wij] is near x = yi. Actually, we 
shall push our calculations up to the level of G^^^ , which 
depends on the first covariant derivative of the curvature 
tensor, and we shall therefore need all the brick potentials 
to be at least G^ near x = yi. 

The application of the above results yields the follow- 
ing explicit expressions for the 1-regular part of the two 
separate 0{mim2) delicate contributions to Wij (defined 
in Eqs. ^^J8^ ~^^ ^ ) 
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(3.32) 
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(3.33) 



As the potential Ri has a source of the same type as 
Wij (namely S{x ~ y^) terms plus a non-compact term 
quadratic in the V potentials) , the calculation of its reg- 
(3.31) ular part can be done in exactly the same way as Wij . Ri 
contains terms oc m^,m2 and mim2. The 0{mi) piece 
is purely singular, the 0(m|) one is purely regular, while 
the 0(7711 7712) one is a mix of regular and singular terms. 
As above, we can decompose the mi 7712 part of Ri in two 
pieces, say 



p[mim2] _ A[mim2] . p[)nim2] 



(3.34) 
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Applying the above results then yields the following 
expressions for the 1-regular parts of these quantities: 



rp [mim2]i 



'12 



(?^l2^'l) 



2(712^1) 



12 



{ni2Vi) ( J_ 
5* 



3(^2^2) 1 / J_ 

2 J V^Vfl 2ri2 

2K2.i)-%^ + 2(n2.0 



3(n2t'2)^ i ( 1 
1*1 



2 r^V^i^^ 



+ 



2ri2 V-S" 



(3.37) 



[k 



[)nim2]i 
(0) i 



'12 



f ("12^2) 

I 2 

ri2 



1 

- 2(rit;2) 



/ 1 



-2(ni2W2) 



3(7112-^1) 



?'2fl2 



2r 



12 



1 



7-2fl2 

(3.38) 



One should substitute the expansions p.28p - p.3ip into 
these resuhs to get their exphcit forms. 

Finahy, we come to the most comphcated 2PN brick 
potential, namely X. It contains contributions propor- 
tional to ml mi m2, mi; m?, m\ m2, mi ml and ml (see 



Eq. (I5TT|) ). The terms in ml, 



'-2J ' 



, m2 are easily dealt 



with (they are either purely singular or purely regular). 
Many, but not all, of the mi m2 terms can be dealt with 
in the same way as the mi m2 terms in Wij and i?^. If 
we again decompose X^™-^™-'^^ in two pieces 



^[mim^] ^ + ij;;]^'"^! , (3.39) 



we have the following results for their regular parts: 



[mim2]i _ 2 

(0) \r - "1 



'^2 



ViSJn ^12 

-(r.i2«if ((^)^- 

3(ni2z;i)(ni2«2) / / 
2 U 

-5(ni2i;2)(riwi) (j^^ - (riwi)^ 
f 2(rit;i)(ri7;2) + ^) ^ ~ (""i"^^' 



1 



1 

7'12 



a). 



-'"""'(K^Ts), 



2ri2 



+ 



ri2 {.Sjjjy 
2(ni2.i)(ri.i)(-^)^ 



'1 

+2(riW2)(n2Wi) (^^r;^^ " ^(riWi)(n2W2) (^^:^ 

[ijor"'^.. = «2^ + - (^) + - ) 

^ ' \r2r12 r2 \S J ri2 \S J 

-{V1V2) ( — 
\r2 



R 



R 

2^, 



7'12 

R 



1 



?'2?'12 



'i{n\2Vi){nx2V\) 



1 

7'2 



-2(nl2^'2)(^^2^'2) (^^^ + 5(ni2Ui)(n2i;2) 



1 

ri2 



1 

?'12 



\^ RJ 



-{n2V2f' 



\^ / R. \^ / R 



1 

f2 



1 

^2 



(3.40) 



f2(7l2W2)(n2t'l) (^(^^ 

t-2(7ii2fi)(rit;2) (^^^^ + 2(n27;i)(riz;2) (;:;^) ^ - ^('^2«2)(riVi) (j^) ^ ■ (3-41) 
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One should substitute the expansions (I3.28p - p.3ip 
into the corresponding terms in Eqs. (I3.40p . (|3.41l) . How- 
ever, these equations involve new types of terms, not dis- 
cussed above. These new terms are of the form 



1 



1 



fi^) ( „2 ^2 + ^3 5 



(3.42) 



-2 0-2 



The fact that we have a specific combination of S 
and rj~'^S'"^ simplify things. Indeed, using the expan- 
sions p.25|) . (|3.26p above we have 



fix) _ fix) 



St 



Sq _r, So 

1 

rf ri 



^0 



02 



(3.43) 



fix) _ fix) 



o4 
•^0 



C3 C2 



5*0 



rl n So S^, + ■ 

(3.44) 

When summing these two equations we see that the terms 
(X l/rf cancel. We shall deal later with these terms, 
which turn out to be delicate to handle but, anyway, 
in the sum of (I3.43P and (13.441) they cancel out. The 
remaining terms either contain an odd power of ri (and 
are therefore purely singular, Eq. (13.171) ) or a positive, 
even power of ri (which makes them purely regular). As 
a consequence, the regular part of the combination p.42p 
reads 




(3.45) 



Note that, thanks to the cancellation of the l/r2 
terms, we have again a property of commutativity 
R[fix) (p{x)] = fix) R[(pix)], for the special type of 
terms ^ix) entering Eq. (j3.42p . 

Concerning the mi m| contribution to X, it is the sum 

of 



^(0) 



and 



X 



[mim2] 
(0) 



1 




1 


(3.46) 


2^?2 


2rirf2 


2rir22 


1 1 


1 


rl 




2r3 ^ 16r'J 
^ rf 


16r|ri 

„2 

'2 


0»,2 3 
^'2' 12 

3 


15ri 


rl 


32r?r22 
ri 


\Qr1rl2 

7-2 
'12 


32r22r22 
(3.47) 


2rlrl^ 


2r|ri2 


32r2r3 ' 



Using the rule (|3.17l) we easily see that each term is 
clearly either purely regular or purely singular. Com- 
puting the regular part of is then easy. 



The most delicate contribution to X is its Oiml 1112) 
one, which can again be written as the sum of 



X 



[771^)712] _ 

(0) 



1 

2r- 



and 



X 



[■m\m2] 
(0) 



1 


1 






I 16r2 


16rlr2 


2rlrl2 






2 
ri 




15r2 


^' 1' 12 


oo„3„2 
OAI 2' 12 


±\ji 2 ' 12 




r-2 
'2 


r2 


r2 
'12 


(3.4 


2rfr22 


2rfri2 


32r2r3 


1 


+ "\ 


1 


(3.4 


^'12 


2r2rf2 


2r2r?2 ■ 



Actually the part x'"^^'-"^^^ is easy to discuss: its reg- 
ular part is obtained simply by discarding the term: 



(0) 



ri/(2 r2 ri^j). Similarly, most of the terms in X 
are easy to treat, being either purely regular, or purely 
singular because of Eq. p.l7p . However, the third, 
fourth, eighth and last terms in the right hand side of 
Eq. (I3.48P are somewhat tricky. [These terms correspond 
to the "dangerous terms" in xC^i^i^^) ^j^g^^ were dis- 
cussed in Section HI when making the link between the 
e-regularization and the Hadamard one.] The third term 
is 



Q = 



1 



16 rf r2 



(3.50) 



while the sum of the fourth, eighth and last terms reads 

15 r2 r?2 



P = -- 



12 



32r2r22 32 rf 



(3.51) 



Both Q and P are of the form f{x)/rl (but we shall 
see that Q is special compared to P). The computation 
of the regular part of f{x)/rl is a bit subtle. It can, 
however, be done by brute force, namely by replacing the 
smooth function fix) by its Taylor expansion around yi: 



fix) = fiyi)+r\dUiyi) + ^,r\r{d,,fivi) 

+ ^, r\ri r?%fe/(yi) + ... (3.52) 

When replacing r\ — > ri n\ and dividing by rl one sees 
that the regular part of f(x)/rl will only come from the 
r'l^*^-'' with £ = 2,4,6,... Moreover, by 



terms r^ 



decomposing rj" — ri in irreducible tensorial parts, as 
in 



11 

r' "* 



■11 
-t ■> 



3 



(3.53) 

= — i (5*-' denotes the symmetric 

trace- free projection of = n\ n\ , we see (in view of the 
definition (|3.15l) ) that only the pieces containing at least 



wher 



n — 
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one Kroneker 6 in the decomposition of nf will contribute 
to the regular part. For instance, in the case £ — 2 only 
the 5^^ in p.52|) will contribute to the regular part of 
f{x)/r\. More generally, we have that R[r^ /r\] ~ [r^ — 
r{)lr\. 

Applying this method yields the following result (here 
written with the simplified notation used around Eq. 
p.lSp ) for the regular part of f{x)/rf, 



6 



A/(0) + — x'9,A/(0) + — A/(0) 



10 



28 



-^rVa.A2/(0) + O(x4). (3.54) 



As one sees on Eq. p.54p (and as can be proven to all 
orders), all the terms on the right-hand side of p.54p are 
derivatives of the Laplacian of f{x) (taken at cc = yi). 
As a consequence, in the particular case where A/(a;) — 
0, the regular part of f{x)/r\ is exactly zero. This is 
the case for the term Q in Xl™?'"^]^ Eq_ ((XsUl) . [Let us 
point out in passing that the discussion in Section III 
above of the link between the £-regularization and the 
Hadamard one essentially consisted in remarking that all 
the "dangerous" terms in xl™!™^! had this innocuous 



structure f{x)/r{ 
have simply 



(2+2£) 



with Af{x) = 0.] Therefore, we 



Qb. = 0. 



(3.55) 



On the other hand, this is not the case for the term P, 
Eq. p.5ip . The evaluation of the regular part of P needs 
to appeal to the result (|3.54p and yields (modulo terms of 
order 0{rf) that will not be needed in our calculations) 



8rh 



' 12 



(rini 



15 
32^^^ 



224 ' 
(rinia)^ 



rl 
15 

3 



' 12 

32 



1 



rl7^l2) 



1'2/ _R 



(3.56) 



Summarizing: We have explicitly displayed all the 
rules needed to compute (near particle 1) the regular 
parts of the various brick potentials V, Vi,Wij, Ri, X en- 
tering the 2PN metric. By replacing V Vr, . . . ,X — > 
Xr, in Eq. p.2p . we define a regularized version of the 
2PN metric generated by two point masses, gj^'^{x) = 

5/ji/[l^fl(a^), . . . ,Xr{x)], which is smooth near particle 1. 



IV. COMPUTATION OF THE INVARIANTS 
ENTERING THE TIDAL ACTION 

As we explained above, when neglecting terms 
quadratic in the tidal parameters /i^^-* etc., the tidal part 



of the two-body action is simply obtained by evaluating 
the 'S'nonminimah Eq. (|2.12p . as a fuuction of the world- 
lines, by replacing the metric gfiu{x) entering the right- 
hand side of (|2.12l) by the (regular part of the) point- 

This reduced 



mass metric gP°™* (x , yi , 2/2 , "t, i , m2 ) 



action is a sum over the various tidal parameters, A*^', 
We can therefore compute separately the 



part of the reduced action associated with each of them. 
This is what we shall do in this Section for the actions 



and 



associated with the parameters ^I'^J^, /^a=i^ '^a=? 
fJ-'\^Si \ [We shall only explicitly write down the results 
for A = 1 but, evidently, they also yield the results for 
A = 2 by exchanging 1^2.] 

First, let us note that each action, say associated with 
the parameter /zi related to the first worldline, is of the 
form 



J dtLf,^{yi,y2,Vi,V2) 



(4.1) 



where the Lagrangian L^^ is the product of a geometrical 
invariant by dri/dt. For instance 

We shall separately evaluate each geometrical invari- 
ant, G^^, G^bc' • ■ before multiplying it by the (regular- 
ized) proper-time redshift factor dri/dt ("Einstein time 
dilation"). Note also that we systematically work with 
the order-reduced 2PN metric, i.e. in which the higher 
time derivatives of yi and 1/2 have been expressed in 
terms of positions and velocities only (yi, 1/2, 'I'l, i'2) by 
iterative use of the (harmonic-gauge) equations of mo- 
tion. As has been discussed long ago, such an order re- 
duction of the action is allowed, when it is understood 
that it corresponds to a certain additional change of co- 
ordinate gauge [45l - l47j . As we shall ultimately be inter- 
ested in computing gauge-invariant quantities associated 
with the EOB reformulation of the dynamics, we do not 
need to keep track of this coordinate change. 

Let us start by discussing the simplest (and physically 
most important) geometric invariant, namely the one as- 
sociated with the electric-type quadrupolar tide, say 

Ja = [Gab Gah]l = [Gc^ G"^]i 

= [Ro.f^D.R^^sxg'^^ g^^u^'u•'u''u\, (4.3) 

where = dy^/dri = (c,Vi) dt/dri. Using two indepen- 
dently written codes (one based on the Maple system, 
and the other one based on the Mathematica one sup- 
plemented by the package xAct we have computed 
the right-hand side of (|4.3I) within the (regularized) 2PN 
metric. [Actually, the Mathematica code regularized, a 
la Hadamard, the value of Ja computed with the full 
(non-regularized) 2PN metric] 

As the PN expansion of the quadrupolar tidal tensor 
(P?7)) starts as 



Gab 



C- Ral 



ObO 
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= + 2 500 - 9ao 360 - dbo gao + doo gab) + ■■ ■ 

one sees that the 2PN-accurate metric (i.e. the knowl- 
edge of goo up to 0{l/c^) terms inchided, of goa up to 
0(l/c^) and of gab up to 0(l/c^)) is exactly what is 
needed to be able to compute Gab to the 2PN (frac- 
tional) accuracy, i.e. Gab = + dabV + c^^{. . .) + c^'*(. . .). 
The same is true for the higher electric tidal moments 
Gabci ■ ■ ■■ However, one can easily see that one loses a 
PN order when evaluating either the magnetic tidal mo- 
ments Hab,Habc, ■ ■ or the time-differentiated electric 

one Gab, The result we obtained, for general orbits, 

is 



Ja 



' 12 



G_ 

ri2 



(5mi 



- 3(ni2i;2) 
6m2) 



3wf 



3(ni2Wi2)^ + 12(ni2W2)^(n-i2Wi2)^ + 6{ni2V2) 

-9t;i2('^12Wl2)^ - 6(ni2Wl2)^(w2Wl2) 
-6(ni2W2)(ni2Wl2)(t'2t'l2) - 3W2("12W12)^ 
-9t;^2('^12W2)^ - 3W2("-12W2)^ + 6^12 
-|-6i;i2("2Wl2) + 3(t;2Wl2)^ + 3W2«12 

Gmi / 109 



("-12W12 



+ Y(f^l2W2) 



ri2 V 4 ^ 

2 



21 



''12 



/365TOf 125TOim2 
TfTV 28 2 



21m^ 



and 



(4.4) 



Similarly, we computed the geometrical invariants 

Jd = [Gabc] 1 = [Ga^-y G"'^^] ^ , 



G(4mi + 10m2) 



''12 

n2/„ „. \2 



lO{ni2V2Y {ni2Vi2)'' + 10(ni2W2)'* 
14 



■y^l2('il2Wl2)^ 



4(ni2Wi2) (W2W12) 



-12l;^2(f^l2W2)^ - 4(71121^2) (^^12^^12)(W2W12) 
-2^2(7112^12)^ - 4w|(r7i2W2)^ 

+ -^^12 + 6Wi2(w2l'12) + 3(W2W12)^ + 3w|wi2 



Grrii 



ri2 



- 32(ni2Wi2)^ + 2(ni2W2)(ni2Wi2) 

16 \ 

+22(ni2«2)' + yt'?2) 

(l2(ni2Wi2)^ + 45(ni2i;2)^ - 18v 



r\2 

n 259miTO2 
— (9"^i + 5 



12 



54r7i; 



(4.7) 



18G2 



' 12 



- («12Wl2)^ +W?2 



("12^12)'^ + 4(ni2W2)^(»^12Wl2)^ 

-3Ui2(»^12Wl2)^ - 2(ni2Wl2)^(w2Wl2) - v\(ni2Vi2f 
-2(ni2U2)(ni2Wl2)(w2Wl2) - 3Wi2("-12W2)^ 



^12 V 
^2 I „,2„,2 



-|-2Ui2 + 2Wi2(w2Wl2) + (W2W12) +"^2^12 
^2G/TOi ^ _ W/„.^„,„\2 „.2 

''12 V 3 



+ nT-2^ (^(rM2Wi2)^ - 



(4.8) 



(4.5) 



The result (|4.4I) . after multiplication by the redshift 
factor 



and 



— - -500 - 2.90.- 



1/2 



(4.9) 



^auRv ^^ks\ g g U U U U 



(4.6) 



which evaluates to (we use again the notation e = 1/c, 
and henceforth often set Newton's constant to one) 



Note the factor i, introduced in the definition of Jb 
to have Jb = (c-R*^^^u'^u'')^, analogously to Ja = 
{RanPuU^u")^ ■ Let us also note in passing that, in 
evaluating Jd, i.e. the square of the electric oc- 
tupole Gap-y, Eq. (12. 8p . it is important to use the 
orthogonally-projected covariant derivative Ifj 
stead of {Gap-yY, one evaluates {Ga/^j)^ where Ga/Sj = 
Sym^,^..^ Vq Rpfiyu , one finds a result which differs 
from (Gap-y)'^ by a term proportional to Jg — {GabY (see 
Eq. (P3T|) below). 

The results for these invariants (along general orbits) 
are 



Jd = 



1 



- 2(7112^12)^ ~ 4(7112^2)^ + 3f 



12 



1 



m2 
''12 



4 ™2 / I . ^2 3 2 

Wi H I -(ni2W2) - -jVl 



ri2 



4{viV2) - 2vl) + -^^{imi + 7712) 
^''12 



(4.10) 



provides the 0{^'"i~^'') piece ( "gravitoelectric tidal 
quadrupole" ) of the reduced two-body action at the 2PN 
approximation level, i.e. including tidal correction terms 
that are (v/c)'^ smaller than the leading order tidal La- 
grangian which is simply given by Ja^"^ — Qm\/r\2- Sim- 
ilarly, multiplying the results (|4.7I) . (|4.8p by the redshift 
factor (|4.10p provides the reduced tidal actions associated 
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with Jd = [G^fccli ^'^d = il-f^abli' the 2PN approxi- 
mation for the electric-octupole terra J^, and at the IPN 
approximation for the magnetic-quadrupole term Jf,. 



In view of their complexity, the resuhs (|4.4p . (|4.7I) . 
(|4.8p . which provide the action for general orbits, are not 
very useful as they are. In the following, we shall extract 
the physically most useful information they contain by: 
(i) focussing our attention on circular orbits, and (ii) re- 
formulating our results in terms of the EOB description 
of binary systems. Note in passing that though circular 
orbits are only special solutions of binary dynamics, they 
are the ones of prime physical importance in many situa- 
tions, most notably radiation-reaction driven inspiralling 
binary systems. 

In the following, we shall therefore restrict our atten- 
tion to circular motions. [However, we shall show below 
how this restricted result can crucially inform the EOB 
description of tidally interacting binary systems.] We 
shall also focus on the relative dynamics in the center of 
mass frame. As we see on Eqs. (lOl) . (|i?7)) . (|i?TUl) . 
the various Lagrangians depend only on the relative po- 
sition yi2 = 2/1 — 1/2 and start depending on (individual) 
velocities only at IPN (for general orbits), and even at 
2PN for the invariants themselves (in the case of circu- 
lar orbits). This implies that we shall not really need to 
use to its full 2PN accuracy the relation between center- 
of-mass variables yp^, y™, ^2^, and relative ones 
2/12, ■^12, namely (in the circular case) 



Vi 



y2 



X2+3 



M 
ri2 c? 



Vl2: 



-Xl+3 I VXX2 

\r\2 



fi2. (4.11) 



and the corresponding velocity relations obtained by 
time-differentiating them, using the fact that y\2 ~ 
T\2 1^12 where ri2 is constant and ni2 rotates with an 
angular velocity given by 



M 
F 



Xi 



1 + £2(^-3) — 
'12 L Tl2 



41 ,\ f M 
4 J \ri2 



Here and below we use the notation 
mi m2 



(4.12) 



v = Xi X2 , X12 — X1—X2 , 
(4.13) 

(recall that M = mi + 7712 so that Xi + X2 = 1). In 
our calculations below the = 1/c* contributions in 
Eqs. (|4.11|) . ()4.12|) do not matter, and can be neglected 
from the start. 

Using such an additional circular (and center-of-mass) 
reduction, we get a much simplified result for the electric- 
quadrupole invariant Jq, Eq. (14.31) . namely 



j(circ) 



6M^X^ 



2 {Xi - 3)M 



ri2 



''^-z^{7l3Xf - 805Xi ~ 336) 



(4.14) 



In a similar manner, one gets much simplified results 
for the other (subleading) geometrical invariants of tidal 
significance, namely the magnetic quadrupolar term J^, 
Eq. (|4.6p . the electric octupolar term Jd, Eq. (14. 5p . and 
also for the time-differentiated electric-quadrupole cou- 
pling. 



^2 



^[iu^V^G^p)iu''V,G"^)]^. (4.15) 



Among these invariants, the 2PN accurate metric allows 
(like for G^^) to calculate to 2PN fractional accuracy only 
the electric-octupole term Jd- The other ones can only 
be computed at IPN fractional accuracy because of their 
"magnetic", or "(?o = c^^dt" character. Our explicit "cir- 
cular" results for Jh, Jd and Je are 



f(circ) 



ISX^M^ 



e' — iiXl 



3ri2 



Xi-9) 



(4.16) 



' 12 



1 + 6^(6X1-7) — 
ri2 



,4^,.l ^2 

irl2 



[QlXi + 4X1 - 98) 



1 + e\Xl - 7) 



M 
ri2 



(4.17) 



(4.18) 



To complete the above results, and allow one to com- 
pute the corresponding associated Lagrangians, let us 
note that the circular value of the redshift factor is 

rfri ^ ^ ^ ^ _ M(Xi-l)(Xi-3) ^2 
dt ~ Fi 2ri2 



AP{Xi - 1) 

Hr^ 
°' 12 



{3X^ - 9X^ + nXi - S)e'^ . (4.19) 



Let us also quote the value of the inverse redshift factor, 
Fi (analog to a Lorentz 7-factor 7 = l/-\/l — v^/c"^), 
namely 



F - - 1 
^ dri 



Af(Xi-l)(Xi-3) ^ 
2ri2 



8^ 



{Xi - l){Xf + 5Xf - 17X1 + 15)e^ . (4.20) 



V. EOB DESCRIPTION OF THE TIDAL 
ACTION 



We have computed above the effective actions associ- 
ated with the tidal parameters ^^^\ ^^^\ a''^^'^ and Hi^'^ ■ 
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Before the restriction to circular motions (in the center- 
of-mass frame) they have the general form 

m J dtL^^{yi2,Vi,V2) , (5.1) 

where fii denotes a generic tidal parameter, and yi2 = 
yi — V2- In this Section we discuss how one can de- 
scribe the actions (|5.1|) within the EOB formalism. Let 
us recall that the EOB formalism jT6l - [l9| replaces the 
(possibly higher-order) Lagrangian dynamics of two par- 
ticles by the Hamiltonian dynamics of an "effective par- 
ticle" embedded within some "effective external poten- 
tials" . For non-spinning '63] bodies the original (velocity- 
dependent) two-body interactions become reformulated 
(and simplified by means of a suitable contact transfor- 
mation in phase space) in terms of three "EOB poten- 
tials": ^(r-off), i?(reff) _and Q{rcf{,P°^)- The first two 
potentials, Aircs) and _B(roff), parametrize an "effective 
metric" 

dsls = g^^(a;off)d<ff (i<ff 

+rls{dOls + sin^ O.s dipls) , (5.2) 

and its associated Hamilton- Jacobi equation, while the 
third potential Q{rcS,P°^) (which necessarily appears 
at 3PN [l^) describes additional contributions to the 
(Hamilton- Jacobi) mass-shell condition 

= + g^^^Mpfpf + Q(rcff,P°*') , (5.3) 

(where ^ = mim2/M = uM is the reduced mass of the 
binary system) that are higher than quadratic in the ef- 
fective momentum jp"^. Following the EOB-simplifying 
philosophy of Ref. [l^ , we shall assume that the third po- 
tential has been reduced (by a suitable canonical trans- 
formation) to a form where it vanishes with the radial 
momentum p'^ . 

In addition, EOB theory introduces a dictionary be- 
tween the original dynamical variables (positions, mo- 
menta, angular momentum, energy) and the effective 
ones. A crucial entry of this dictionary is a non-trivial 
transformation between the original "real" energy, i.e. 
the value of the original (total) Hamiltonian H , and the 
"effective" energy —p'^ = Hcs entering the mass-shell 
condition (15. Sp . Because of this transformation, the final 
EOB-form of the (original, real) Hamiltonian reads (here 
we set c — 1 for simplicity) 



H''''''{x,s,Pc«) = M^l + 2,.(^^-iy (5.4) 
where HcS = i?efi(iCeff,Poff) is given by 



(5.5) 



Here JcS = XcS'xPcS denotes the effective orbital angular 
momentum, which, by the EOB dictionary, is actually 
identified with the original, total (center of mass) orbital 
angular momentum J of the binary system: Jeff = J. 

Let us now discuss what are the various possi- 
ble methods for reformulating an original action of 
the type Lo{yi2,Vi,V2, ■ . ■) + ^il Lf^^{yi2,Vi,V2) [where 
fii stands for a sum over a collection of tidal pa- 

, (2) (2) (3) (3) 1 . , 

rameters , ^J'2 % /^i % /^2 % • ■ -J i™o correspondmg ^i- 
deformations of EOB potentials: Ao{rcs) + Mi ^/Ji {fcs), 
-Bo(roff) -I- /ii i3^i(rcff), Qo{rcS,p°^) + fJ,iQf,^{rcB,p''^)- 
The main difficulty in finding the perturbed EOB po- 
tentials A^j^ , _B^j and Q^-^ that encode the dynamics of 
is that such a dynamical equivalence is obtained only 
after some a priori unknown phase-space contact trans- 
formation between the EOB phase-space coordinates, say 
^off = (a;off,PofT) and the original (harmonic-coordinate- 
related) ones, say £,h — (2/12, 'Ui2)- For simplicity, we 
assume that we have already performed the reduction of 
the original harmonic-coordinate dynamics to its center- 
of-mass version, in which one can express Vi and V2 
in terms of the relative velocity V12 = Vi — V2 and of 
2/12 = 2/1—2/2- On the other hand, we do not imme- 
diately assume that the original Lagrangian dynamics is 
expressed in Hamiltonian form. [Let us recall that, as 
was found long ago 0, HH , starting at the 2PN level, 
the harmonic-coordinate dynamics does not admit an or- 
dinary Lagrangian L(y,y) but only a higher-order one 
L{y,y,y). In order to express the 2PN dynamics in 
Hamiltonian form one already needs some (higher-order) 
contact transformation. However, this transformation is 
well-known, e.g. [46| . and we do not need to complicate 
our discussion by explicitly mentionning it. It will be, 
however, taken into account in our calculations below.] 

The transformation T between ^ofi and will have the 
general structure 

£,h = To(Ccff ) + Ml ($off ) ■ (5.6) 

The unperturb ed p art T'o(Coff) is known from the pre- 
vious EOB work [ifllll, but the 0(mi) perturbed part 
T)ii('Ccff) is unknown, and, actually, is part of the prob- 
lem which must be solved for reformulating the (per- 
turbed) harmonic-coordinate dynamics in EOB form. 
This means, in particular, that it would not be correct to 
try to compute A^-^ , and , by simply replacing in 
the tidal action (|5.1|) the harmonic variables by their 
unperturbed expression To{^cs) in terms of the effective 
variables ^off- 

For the general case of non-circular orbits, a universal, 
correct method for transforming the original Lagrangian 
L{£^h) = Lniih) + Mi-^Mi(C/i) in EOB form consists (as 
explained in flSl) of the following steps: (i) to trans- 
form the original Lagrangian L{£^h) in Hamiltonian form 
H{S.h) = Ho{^h) +Mi^^Mi(Cff), where £,h = iq,p) are 
canonical coordinates; (ii) to extract the gauge- invariant 
content of H($^h) by expressing it in terms of action vari- 
ables la = § Pa dqa, which yields the Delaunay Hamil- 
tonian H{I) — Hq{I) + /if Hfj_-^{I); (iii) to do the same 
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thing for the EOB Hamiltonian, i.e. to compute, as a 
functional of the unknown EOB potentials, its Delaunay 
form ffEOB(/) = H^'^^il) + fn HJ^^b^j-j. ^^^^^^ 
(iv) to identify the known H{I) to Hbob{I), which de- 
pends on the unknown functions A^-^ , B^-^ , . This 
last step yields (functional) equations for , B^-^ , Q^j^ 
and thereby allows one to determine them. [In prac- 
tice, the functional dependence on A,B,Q is replaced 
by a much simpler parameter-dependence by using the 
method of undetermined coefficients for parametrizing 
general forms of A,B,Q.] An alternative (and equally 
universal) method for transforming L{S,h) in EOB form 
(as used in Q^) is to add the transformation = T{S,cs) 
to the list of unknowns (using the method of undeter- 
mined coefficients), and to directly solve the set of con- 
straints for T, A, B and Q coming from the requirement 
that HEOBiicS, A, B,Q) ^ H(T{^^s)). [One must then 
take into account that T is constrained to be a canonical 
transformation.] 

The 2PN-accurate results, given for several tidal in- 
teractions in the case of general orbits, in the previ- 
ous Section, can, in principle, be transformed within the 
EOB format by using any of the two methods we just ex- 
plained. However, from the point of view of current astro- 
physical applications, one is mainly interested in knowing 
the EOB description of (quasi) circular motions. In this 
case, we a priori know that it is only the A radial po- 
tential which matters. Knowing this, the question arises 
to compute the tidal perturbation A^^-^ of the EOB A 
potential in the most efficient manner; possibly without 
having to go through the rather involved, general univer- 
sal methods recalled above. Fortunately, it is possible to 
do so by using the following facts. 

The first useful fact concerns the relation between the 
tidal perturbation (in harmonic coordinates) of the La- 
grangian of the binary system, say 

5L^(yh, Vh) = Ml i^i {yi2,Vi,V2) (5.7) 

to the corresponding tidal perturbation (in harmonic- 
related phase-space coordinates) of the Hamiltonian, say 

SH''{yh,Ph) = Hty^{yh,Ph) - H^ia^i_f^^^{yh,Ph) ■ (5.8) 

[Strictly speaking, as we recalled above, the harmonic- 
related (g,p) = (jjhjPh) phase-space coordinates involve 
a supplementary 0(l/c'*) gauge transformation linked to 
the order reduction of L2PN{y,y,y) into U^pNiv' 'V')-] 
Note that here and in the following the notation SQ{£^) 
will always refer to the tidal contribution to some func- 
tion of specified variables, i.e. SQ{^) = (5fuii(0 ~ 
'3tidai-frco(0- One has to be careful about which variables 
are fixed as, for instance, the transformation between La- 
grangian {q,q) and Hamiltonian {q,p) coordinates does 
contain a tidal contribution (because S*^^'^^^L{yh,yh) does 
depend on velocities). This being made clear, we have 
the well-known universal result about first-order defor- 
mations of Lagrangians by small parameters, L{q, q) = 



Lo{q, q) + Ml (<?' 9), namely 

dH'\yh,Ph) = -5L''{yh,Vh) (5.9) 

which follows from the properties of the Legendre trans- 
form. 

Let us now apply the second method recalled 
above for transforming the "harmonic" Hamiltonian 
Htii{yh,Ph) = H^{yh,Ph) +5H^{yh,Ph) (where the in- 
dex refers to the unperturbed, tidal-free dynamics) 
into its corresponding EOB form i?fuii^(xEOB,PEOB), de- 
fined in Eqs. (15. 4p . (15. 5p above. [For clarity, we denote 
here the effective-one-body phase-space coordinates by 
s^eoBjPeob, instead of XaS,Pcfi as above.] The crucial 
point is that the EOB potentials entering the definition 
of Hf^^ must be the full, tidally-completed values of 
A, B and Q, e.g. 

^full(rEOB) = ^o(rEOB) +Mi^mi(^EOb) 

= Ao(rEOB)+M(rEOB). (5.10) 

In other words (5i?^°^(a;EOB,PEOB) = 
■^fun^(2^EOB,PEOB) - -ff(?°^(a;EOB,PEOB) is ob- 
tained by varying the functions A, B and 
Q (i.e. ^(rEos) = ^o('"eob) + <5^('^eob), 
etc.) in the definition ((5^ . (I53|) of 

Hf^i^ [X^OB, PEOB ; ^(r-EOB ) , ^ (rEOB ) , Q (^EOB , )] ■ 

This second method for mapping -f^fun(?/i) 
into Hfun^(^EOB) [where = {yh,Ph), ^eob = 
(s^eoBiPeob)] consists in looking for a full, i.e. per- 
turbed, (time-independent) contact transformation 
C/i — 7tuii(CEOB) = 7o(^eob) + Mi^a'i(^eob) that 
transforms i?fuii(6i) into H^^^ {(^eob) , i-e. such that 

HUTfuni^EOB)] = i/^°^(&OB) . (5.11) 

Rewriting the full transformation Tfun as the composi- 
tion T' o Tq of the known unperturbed (tidal-free) con- 
tact transformation ^JJ = To(i^eob) mapping (^°) i^to 
Hq'~'^ {£^eob) with an unknown near- identity additional 
transformation, ^ = T'(^O)^ = + Mi{G^,(C^),?^} 
[where {/, 5} denotes a Poisson bracket and where 
MiG^i(Ch) i^ the first-order generating function associ- 
ated with the canonical transformation T'] , and expand- 
ing all functions in Eq. (|5.1ip into unperturbed plus tidal 
contributions {H'' ^ H^} + 5H^ , T = (1 + 5T') o Tq, 
jjEOB ^ jjEOB ^ ^^eoB) ig^jjg ^i^g condition 

= <5i/^=OB(eEOB), (5.12) 

where 5G(^0) = miG^, (^D- 

In general i5G(^^) is part of the unknown functions that 
must be looked for when writing the condition (|5.12p . 
However, another simplifying fact occurs in the case 
where one focusses on circular motions: the supplemen- 
tary term {SG, H^} happens to vanish. Indeed, SG{£,1) 
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is a scalar function and the Poisson bracket {6G,H'^} 
is equal to the time derivative of SG{£,'^) along the H'^- 
dynamical flow, which clearly vanishes along circular mo- 
tions. This allows one to conclude that, along circular 
motions, we have the simple condition 



C?.=To(eEOB) 



6H 



EOB 



(^eob)]' 



(5.13) 



where the left-hand-side is, in principle, fully known. 

Let us now evaluate the right-hand side of Eq. (|5.13l) . 
When restricting the definition ([Q]) . ([53]) of the EOB 
Hamiltonian to circular motions, the terms {pf~'^)'^ / B 
and Q{rEOB,P^'^^) disappear (because one works with 
a gauge-reduced Q which vanishes with pf'~'^). As a 
consequence -ffEOB (''eob , J) only depends on the A po- 
tential. The difference <^-ffEOB = -f^EOB I'^eob , >/, ^fuii] — 
HeobI'^eob, J, Aq] can then be simply computed by vary- 
ing A (Afuii = Aq + SA) within H^'^^[A]. To write explic- 
itly the result of this variation it is convenient to work 
with dimensionless variables. We can replace the two 



phase-space variables Teob, P^*^^ = J that enter i^EOB 
by their dimensionless counterparts 

_ GM _ G(mi + TO2) 



TEOB c2 TEOB 



and 



c J 



cJ 



GM fi Gmim,2 



(5.14) 



(5.15) 



In terms of these variables the explicit expression of 
reads 



= Af ^1 + 2 1/ (-1 + x^A{u){l + f . (5.16) 

Varying A{u) in (|5.16p then yields the following explicit 
expression for the right-hand side of (|5.13p : 



J 



1 + p 



^ A{u) [1 + 2 1/ (- 1 + ^A{u){l+fu^)'^ 



SA{u) 



(5.17) 



In addition, one must take into account the constraint 
coming from the reduction to circular motions, namely, 
from = -dH^^^/drEOB, the fact that du[A{u){l + 
P V?)] = 0, i.e. the fact that p is the following function 
of u (using a prime to denote the u-derivative) 



„-2 _ „-2 

J 7circ 



{u) 



A'{u) 
{u^A{u))' 



(5.18) 



Note that this relation depends on the value of the radial 
potential A{u). If one is considering the full, tidally- 
perturbed circular motions one must use ^fun(w) = 
Aq + 5 A in Eq. (I5.18|) . On the other hand, as we are now 
interested in considering the (first-order) tidal perturba- 
tions 6H'^ and 5H^°^, and their hnk in Eq. ((Sl^ . we 



can evaluate SH^^^^ with sufficient accuracy by replac- 
ing in the coefficient of SA{u), on the right-hand side of 
Eq. (|5.17p . A{u) and by their unperturbed, tidal-free 
expressions Aq{u) and j^g(u) (i.e. by replacing A ^ Aq 
on the right-hand side of (15.181) '). 

Combining our results (|5.9p . (|5.13p and (I5.17P we fi- 
nally get a very simple link between the tidal variation of 
the harmonic-coordinate Lagrangian SL{yh,Vh) and the 
corresponding tidal variation 5A{u) of the EOB A poten- 
tial, namely 



SA{u) = - 



where 



M V & 



VW) [SL{yu,vu)f:ilT,(^^) , (5.19) 



A{u) 
1 + p u? 



l+2v{-l 



+ ^ A{u){l + p u^)) 



(5.20) 



A=Ao 



Here, the superscript "circ" means that p must be re- 
placed by jcirc('")' Eq. (|5.18l) . In addition, as said above, 
one can use as value of A{u) its unperturbed (tidal-free) 
value Aq{u), both in the explicit occurrences of A in 



(|O0l) . and in the definition ([5?T8| of j^ircl")- Finally, 
if we introduce the short-hand notation 



1 



A{u) = A{u) + -uA'{u), 



(5.21) 
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F{u), Eq. (j5.20l) . can be rewritten in the exphcit form transformation 



,ADM 
12 



^EOB 



one could think of using 



F{u) = A{u) 



l + 2v 



A{u) 



A{u) 



(5.22) 



As we said, it is enough to replace A{u) — A^){u) in the 
above results. In addition, as we computed 8L only to the 
2PN fractional accuracy, it is sufSciently to use a value 
of F{u) which is also only fractionally 2PN-accurate. 
One might a priori think that this would mean using 
for A{u) — j4o(w) in (|5.22l) the 2PN-accurate expres- 
sion of the point-mass A potential, namely A^^{u) — 
1 — 2u + 2vv?. However, the contribution 2vu^ = 
2v{GM/ {c? r-E^o^))^ is 0(l/c^) compared to one, which 
is the leading-order value of F{u), that starts as F{u) — 
1 + 0{u) — 1 + 0(l/c^). The same consideration ap- 
plies to Aiu). [The situation would have been different 
if F{u) had been, say, oc A'{u)] This means that, at the 
2PN fractional accuracy, we can use the value of F{u) ob- 
tained from the leading-order, "Schwarzschild-like" value 
of Ao(w), namely A}^^{u) — l — 2u. The corresponding 
A function is then: Aq^^{u) = 1 — 3 u, so that 



F 



2PN 



(m) = (1-3w) 



H-2i^ -1 



1 - 2ii 



Vi~~3w. 

(5.23) 

Consistently with the fractional 2PN accuracy, and re- 
membering, that u = 0(l/c^), we could as well use the 
2PN- accurate series expansion of (I5.23p . say F^^^(u) — 
1 + fi{>^) " + /2(j^) u"^ + 0{u^). However, it is better to 
retain the information contained in (j5.23p that, in the 
test- mass limit ^ (where Aq (u) — >■ 1 — 2 u) the exact 
value of F{u) becomes 1 — 3u (see later). 

There remains only one missing piece of information 
to be able to use our result (|5.19p for computing the 
various tidal contributions to A{u). We need to work out 
the explicit form of the unperturbed transformation Tq 
between teob and r^. 

A first method for getting the transformation Tq (at 
2PN) is to compose the transformation — )• ^adm (ob- 
tained at 2PN in and at 3PN in ^) with the trans- 
formation iJadm — >■ ^EOB (obtained at 2PN in [l^, and 
at 3PN in ^8j]). For our present purpose, it is enough to 
restrict these transformations to the circular case, i.e. to 
transformations rh — >■ tadm and tadm — ^ r^oB- 

The transformation h — ^ ADM starts at 2PN, 

i.e. y'X = a^^DM ^ ^-4y2PN(3.ADM^pADM)^_^ ^i^l^ 
y2PN(3.ADM pADM) gj^g^^^ g_g^ gq. (4.5) of jS^]. ItS 

circular, and center-of-mass, reduction (with ni2-pA = 0, 
Pi = P2^P and (p//i)2 = GAf/r 12 + 0{l/c^)) yields 
at 2PN 



' 12 



^ADM 
^^12 



29 



GM 



(5.24) 



On the other hand the transformation ADM — >■ EOB 
starts at IPN. To determine the corresponding radial 



Eq. (6.22) of [lOi]. However, this equation needs to be 
completed by the knowledge of the circularity condition 
relating (p^^^/fi)'^ to GM/r^^^ at the IPN level in- 
cluded. This IPN-accurate circularity condition can, e.g., 
be obtained from combining the IPN-accurate r^^^ — 
^ADMf^j^ relation given in [5i[ (see below), with the fact 
that (setting uadm = GM/(c^ tadm)) (Padm/(mc))^ = 
J^^ADM- This yields {pADu/ilJ-c))'^ = wadm +4u^o]^, 
and therefrom the relation between tadm and rgoB • An- 
other method (that we have checked to give the same 
result) for determining the rf^^ — >■ teob transforma- 
tion, which does not use Eq. (6.22) of [ig], consists in 
directly eliminating the dimensionless angular momen- 



tum j between the two relations 



„ADM 



„ADM 



,EOB 



U) and 



^'-'^(y). The former relation was derived at 



3PN in Ref. [511 and reads, at 2PN, 



„ADM 



GM 



■J 



4 1 1 ■ 

l----(74-43^)- 



(5.25) 



while the latter one is obtained by inverting the 2PN- 
accurate version of Eq. (I5.18p . namely, using v42pn(u) = 
l-2u + 2i^u^: 



Inserting ([g?^ into ((05|) yields (at 2PN) 



GM 

„2 „ADM 
^ '12 



1 
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(5.26) 



(5.27) 



Then, combining Eq. ([07)) and Eq. ((OI)) yields the 
looked for transformation r^*-*^ — > ri2 , at 2PN accuracy. 



GM 



' 12 



„EOB 



l + 6iy 



GM 



(5.28) 



or, setting Uh = GM/{c^r\2) by analogy with 
GAf/(c2 teob), 



Uh 



1 



(5.29) 



We have written the transformation (|5.28p . (|5.29p so as 
to exhibit the exact form of the transformation r/j — > 
teob in the extreme mass ratio limit — > 0, namely 
rh = fEOB - GM/(? + 0{v). 

Summarizing: The (first-order) tidal contribution 
5A{u) — iiiA^^{u) to the main EOB radial po- 
tential, associated with any tidal parameter {— 
\ . . .), is given in terms of the correspond- 
ing harmonic-coordinate tidal contribution to the ac- 
tion £Z^/i, w^) = fiiLf^^{yh,Vh), for circular moti on, b y 
Eq. dog, where F{u) is given (at 2PN) by Eq. ([STM)) . 
and where the transformation between the harmonic 
radial separation ri2 and the EOB radial coordinate 
TEOB = GM/{c^ u) is given by Eqs. or ([05)) . 
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VI. EOB DESCRIPTION OF TIDAL ACTIONS 

We have explained in the previous Section how to con- 
vert each contribution ~ fii L^^{yh,Vh) to the (reduced) 
tidal action into a corresponding additional contribution 
jiiA^^iu) to the main EOB radial potential A{u). For 
instance, if we consider the dominant tidal parameter, 



i.e. the electric quadrupolar one, (or ^i2~^' , after 

exchanging 1 2), the combination of the result 
for the associated Lagrangian, with Eq. (|5.19p yields 



,(^=2) 



(6.1) 

In other words, apart from a (negative) numerical coeffi- 

cient, and the rescaled tidal parameter ' /{Mv) (where 
Mv = fj, = mi rail (jn\ + mi) is the reduced mass of the 
system), the corresponding tidal contribution to A(u) is 
the product of three factors: ^ F(u), dri /dt and the ge- 
ometrical invariant associated with the considered tidal 
parameter, e.g. [Ga^G"'']i for the electric quadrupole 
along the first worldline. In addition, two of these fac- 
tors, dri/dt and the geometrical invariant must be re- 
expressed as functions of the EOB coordinates by using 
Eq. Km . 

Let us start by applying this procedure to the domi- 
nant tidal action term: the electric-quadrupole one (|6.ip . 
We have given above, in Eq. (|4.14l) . the 2PN-accurate 
value of Ja = [Gai3 G"^]! in harmonic coordinates. Us- 
ing the transformation (|5.29p to replace l/rj'j in terms 
of 1 /rEOB leads to 



Kcirc) _ 6M^X| 



1 + e- 



[Xi + 3)M 



+e 



r^OB L ''EOB 
.4 ,on.^2 



^^''eob 



{295X1 - 7Xi + 336) 



.(6.2) 



In addition the reexpression of the time-dilation factor 
dTi/dt, Eq. (|4.19l) . in terms of l/rsoB yields 



dri 



1 



1 



-j^ = — = l--iX,-l){X,~3)i 



dt 



lu\Xi - l){Xl - 3X^ + 3Xi + 3)e^ . (6.3) 
8 



Their product yields the electric-quadrupole tidal La- 
grangian (stripped of its prefactor j /i^' ) in EOB co- 
ordinates, at the 2PN accuracy, namely 



G: 



M4 



r 



(6.4) 



Adding the further factor ^y F{u), as well as the prefac- 
tor, leads to the corresponding contribution to the EOB 
A potential, namely 



(2) /I / \ .l(2)LO / N J(2) / N 

fi\ 'Ajr,{u) = A\>^^^^^.^{rEOB)A\>^^^^^.Ju) , 



where 



.(2)LO / X 
dielectric (''EOB ) 



3G2 /if^Af 



'eob 



(6.6) 



(6.7) 



and 



^feLtric(") = Vn^)^a = 1 + a> + ^ 0{u^) , 

(6.8) 

with 



337 o 1 

X? + -Xi + 3 . 

28 ^ 8 



(6.9) 
(6.10) 



The leading-order (i.e. Newtonian level) A potential (|6.7p 
is equivalent to Eqs. prB| . pTTl) above (i.e. Eqs. (23), (25) 
of [a]), using the link 



1 



(2^-1)!! 



(6.11) 



The term of order u (i.e. IPN) in the relativistic ampli- 

^ ('2) 

fication factor eiectric(^)' l|6.8p . coincides with the 
result computed some time ago (see Eq. (38) in 5]). By 

contrast, the (2PN) term of order in ^i^cicctric(") 
the main new result of our present work. Let us discuss 
its properties. 

Similarly to the IPN coefficient ~ ^ Xi which was 
positive, and monotonically increasing (from to 5/2) as 
Xi = mi/M varies between and 1, the 2PN coefficient 
is also positive, and increases as Xi varies between 
and 1. When Xi = (i.e. in the limit mi <C mi) ai 
takes the value -I- 3, while when Xi = 1 (i.e. in the limit 
TOi ^ mi) it takes the value 

849 

a?(Xi = 1) = = 15.16071429 . (6.12) 

56 

Note that this is about five times larger than its value 
when Xi = 0. Of most interest (as neutron stars are 
expected to have rather similar masses ~ 1.4 M0) is the 
equal-mass value of ai, which is 



Xi 



85 
14 



6.071428571. 



(6.13) 



where 



Ca = 1 - -u(X2 - 6X1 - 3)e2 



(6.5) 



56 



(21X1* - lUXf + TUXl + 238Xi + 357)e'' . 



In other words, the distance-dependent amplification 
factor of the electric quadrupole reads, in the equal-mass 
case 



A 



(2) 

1 electric 



ic(«) 



equal-mass 



1 



85 

14' 



0(m3 



22 



= 1 + 1.25 u + 6.071429 zi^ 
+0{u^). (6.14) 

We will comment further below on these results for 
^1 electric ("") ^^^'^ recBut comparisous between nu- 

merical simulations and the EOB description of tidal in- 
teractions. For the time being, let us give the correspond- 
ing results of our analysis for some of the sub- leading tidal 
interactions. 

The EOB-coordinate value of the electric octupole, 



j: 



(circ) 



^ , Eq. dm?]), reads 

(circ) 



7^ 



l-He2(6Xi + l) 



EOB 

•"eob 



M 

fEOB 



{%-iXl + 14Xi + 17) 



(6.15) 



its corresponding action (stripped of its prefactor) is 



from to Xf^™ = 311/1760 = 0.1767046, before increas- 
ing as X\ goes from X™'" to 1, to reach the final value 
Q!^(l) = 211/8 = 26.375 for = 1. Note, however, that 
when (as expected) the two masses are nearly equal the 

factor ^1 ciectric('^) amplification factor. In particu- 
lar, its equal-mass value is 



^lc'lectric(^) 



equal-mass 4 257 o ,0, 

= l + y" + ^"' + 0(w3) 

= 1 + 1.75 m + 5.3541671*2 
+0{u^) (6.23) 



which is similar to its corresponding quadrupolar coun- 
terpart Eq. (|6.14p . 

Let us finally give the corresponding results for the 
magnetic quadrupole and time-differentiated electric 
quadrupole. For the magnetic quadrupole (at the IPN 
fractional accuracy) we found 
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while the corresponding contribution to the EOB A po- 
tential reads 



^fite°tric('^EOB)^i' 
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where 
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and 
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1 electric 



with 



[u) = ^F{u) Cd = l + aiu + aiu^ + 0{u^) , 
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Here, both results ([OT]) and ([6:221) are new. Note that, 
contrary to the quadrupolar case where ai and a2 were 
always both positive (so that A\ ' (u) was always an am- 



plification factor) the electric-octupole factor A^ 



(3) 



:lcctr 



(u) 



is smaller than 1 (for large separations) when Xi < ^ ~ 
0.2667. Moreover, while the Xi-variation of af is mono- 
tonic (going from — 2 to ^■ as Xi increases from to 
1), af(Xi) first decreases from a^{0) = | = 2.666667 
to a^(X™'") = 42853/28160 = 1.521768 as Xi increases 



/:6 = l + -(Xi+3)(3Xi + 5)ue2 
6 



(6.26) 



^fmag„etic(") = Ct = 1 + a\ u + 0{u^) , (6.27) 

with 



ai = Xi 
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Here a\{Xi) is always positive, and monotonically in- 
creases from a\{fd) = 1 to aj(l) = f = 3.833333; its 
equal-mass value being a\ (5) — ^ — 2.166667. 

Finally, for the time-differentiated electric quadrupole, 
we got 
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A^^l{u)^^F{^t,^l + alu + 0{u'), (6.32) 
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with 



a? = -(Xi+2)(2Xi-l) 



(6.33) 



One of the characteristic features of the EOB formal- 
ism for point-mass systems is the natural incorporation 
of the exact test-mass limit v ^ 0. Indeed, in this limit 
the effective metric (|5.2p describing the relative dynam- 
ics reduces to the Schwarzschild metric: lim^^o A{u) — 

1 — 2u = (\imi,^o B{u)^ . Let us study the test-mass 
limit of tidal effects, with the aim of incorporating it sim- 
ilarly in their EOB description. When considering the 
nonminimal worldline action of particle 1, the simplest 
test- mass limit to study is the limit mi/m2 — > 0. [When 
considering tidal effects within body 2, the permutation 
1 O 2 of our results below allow them to describe the 
limit 1112/1711 — > 0. We leave to future work a study of the 
limit 7712 /toi 0, when considering tidal effects taking 
place within body 1.] In the limit investigated here, one is 
considering a tidally deformable test-mass {mi, \ . . .) 
moving around a large mass TO2 ^ wi . The effective ac- 
tion of body 1 is then exactly obtained by evaluating the 
A — 1 contribution of the general (two-body) effective 
action (I2.12p within the background metric generated by 
the (non-tidally deformable) large mass 7712, at rest, i.e. 
within a Schwarzschild metric of mass 7712- The latter 
reads 



Gm,2 
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in "Schwarzschild" , or areal, coordinates, and 



ds'^{m2) = - 
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in harmonic coordinates: ~ — Gm2/c^- As a check 
on the results below (and on our codes), we have com- 
puted them both in Schwarzschild coordinates and in har- 
monic ones. 

The geometrical invariants J a = G^^, etc. take the 
following values in this Schwarzschild limit, and when 
considering as above circular motions (we again set G 
and c to one for simplicity): 
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where us = G 7712/(0^ r^). 

We have indicated in these results the expansions in 
powers of the inverse harmonic radius Vh as checks of 
our 2PN-accurate results, written in harmonic coordi- 
nates, see Eqs. (j4.14p . (|4. 16^ ^ (14. 18p . In the following, we 
shall focus on the transformation of the exact test-mass 
geometrical invariants above into corresponding contri- 
butions to the EOB A potential. As explained before, 
Eqs. (I5TT9)) . (jo]) , apart from the universal prefactor 
— 2/{Mvc^), and the specific original tidal coefficient 
multiplying the considered geometrical invariant (such 

as for the electric quadrupole), the contribution to 

Aiu) associated with some given invariant is obtained by 
multiplying it by two extra factors: (i) the time-dilation 
factor dri / dt and (ii) the EOB-rooted factor F{u). Let 
us discuss their values in the test-mass limit ttji ^ m2 
that we are now considering. 



The first factor is the square-root of 



dri 
~dt 



2Gm2 



dip 



(6.40) 



Denoting, as above, us = Gm2l{(? r^), and us- 
ing the well-known Kepler law for circular orbits in 
Schwarzschild coordinates, ft^ = Gm2/r^, simply yields 



dri 
'dt 



VI - 3 us • 
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The exact test-mass limit of the second factor is ob- 
tained by taking the limit — > in the exact expres- 
sion (|5.22p . In this limit, A{u) — ^ 1 — 2u, so that 
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A{u) 1 — 3w, and 



(6.42) 



In addition, as the EOB coordinates reduce to 
Schwarzschild coordinates in the test-mass hmit — > 0, 
and M ~ mi + 7712 — > m2 , we have simply 



us 



Gm2 



GM 



(6.43) 



^ Ts ^EOB 

In other words, the two extra factors (|6.41|) . (|6.42p be- 
come both equal to y/1 — 3u. As a consequence the 
A contribution corresponding to the various geometri- 
cal invariants (|6.36p - (|6.39l) is obtained (apart from a 
constant prefactor) by multiplying these invariants by 

(Vl — 3 u) = 1 — 3u = 1 — 3^5. Including the uni- 
versal factor —2/{Mi/c^) and the various tidal coeffi- 
cients i J, i ji ... (as well as the factor 4 in 

^ab ~ ^Jb) yields the following exact, test-mass contri- 
butions 
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c2 mi r|oB 



1 -3u, 
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(2) .test-mass. ^ _ 24 — ^ ^'^^ 

'^1 mi r^QB 1 - 3 u 

(6.46) 

'(2) , test-mass. ^ q f^l^^ (™2)^ (1 - 2 m)^ 

Ml ^ '(2) l"; — — y n • 

Ml mi r|QB 1 - 3 u 

(6.47) 

One easily sees that the various exact, test- mass am- 
plification factors A{u) exhibited here, are compatible 
with the Xi — 7> limit of the 2PN-expanded ones ~ 
1 + aiu + (12 + 0{u^) derived above. 

A striking feature of all the amplification factors 
present in Eqs. (|6.44p - (l6.47l) . such as 

9 

J (2) test-mass 



test-mass/ X _ i , o 
^1 electric [U) — i + 6 



1 - 3u 



is that they all formally exhibit a pole cx 1/(1 — 3 u) 
mathematically located at 3u = 1, i.e. corresponding 
to formally letting particle 1 tend to the last unstable 
circular orbit, located at 2>G •012/0^ ("hght ring" orbit). 
This behavior has a simple origin. 

The invariant simplest to consider in order to see this is 
J a = G\j^. From Eq. (|4.3I) its covariant expression reads 



Let us study its mathematical behavior in the formal 
limit where particle 1 tends to the light ring orbit. 
Using the language of Special Relativity, we consider 
the Schwarzschild coordinates as defining a "lab-frame" . 
With respect to this lab-frame, particle 1 becomes ultra- 
relativistic as it approaches the light ring. More pre- 
cisely, near the light ring the lab-frame components of 
the four- velocity = (dt/dTi){c,v^) tend towards infin- 
ity proportionally to dt/dri = Fi = 1/^1 — 3u, while 
the lab-frame components of Ra^ipi, (and the metric) 
stay finite. As G^^ is quartic in the lab-frame com- 
ponents of u^, it will tend towards infinity like T\ = 
{dt/driY — ~ 3u)~^. The corresponding contribu- 
tion to A{u) is obtained by multiplying G^^ by the factor 
{dri/dtY = Fj"^ = (1 — 3m)"'"^ which reduces the blow-up 
of Gl^ to the milder (1-3 u)-'^+^ = (1 - 3 m)-^ behavior 
that is apparent in Eqs. (|6.44p or (|6.48p . 

A different way of phrasing this result uses the law 
of transformation of the electric and magnetic compo- 
nents of the Weyl tensor. Gab and Hab under a boost. 
Using, for instance, the fact that, under a boost with ve- 
locity /3 — tanh ip in the x direction, the complex ten- 
sor Fab — Gab + i Hab Undergoes a complex rotation 
of angle = iip m the yz plane [s^ , one easily finds 
that the transverse traceless components of Fab (in the 
yz plane) acquire, under such a boost, a factor of order 
cos^ -0 — cosh^ (j) — {1 — fP)~'^ = T\. Because of the spe- 
cial structure of the tensor Fab oc diag (— 1, — 1, 2), with 
the third axis z labelling the radial direction, this rea- 
soning shows that boosts in the radial (z) direction leave 
Fab invariant. However, we are mainly interested here 
in boosts in a "tangential" direction, say x, associated 
with the fast motion of a circular orbit, and therefore 
orthogonal to the radial direction, which do introduce a 
factor Fi in some of the boosted components of Fab- For 
completeness, let us indicate that, because of this spe- 
cial structure of Fab, the invariant Ja — G^^ for general, 
non- circular orbits is equal to 



"^a — ^ab ^ 
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1 + 3u,^ 



3ut 
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where u^g = r^((u' 



+ sin^ 0{u'''f) is the square of the 
part of the 4-velocity u'^ that is tangent to the sphere. 
[The radial component of the 4-velocity brings no contri- 
bution to Ja-] 

The behavior near the light ring of the magnetic- 
quadrupole invariant Jf, = jH^f^ is understood in the 
same way as that of Ja — G^^i^- Concerning the other 
invariants, one can note that Jd 
as the sum 



^abc '^^^ be written 



Jd ~ Gafj^ — CaPj G"^''' -I- -— Je 
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where 



(6.49) 
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and 



with 



Gap 



3.53) 



3.54) 



Similarly to G^^, Eq. ([OO)) . the term C^^^ in Eq. (|63T|) 
is quartic in m'' and is therefore expected to blow up like 
r^. On the other hand, though Gai3, Eq. (|6.54p . is cubic 
in it^, it only blows up like (so that Jg ^ Tf and Jd ^ 
C"^ + Je ^ rf ) because of the special geodetic-precession 
properties of the proper-time derivative operator V /dr — 
u^^x (see, e.g., Section 3.6 of [53]). 

Having understood that the formal pole-like behav- 
ior, ~ (1 — 3u)~^, in the test-mass limit of the electric- 
quadrupole A potential is linked to simple boost proper- 
ties of Gab near the light ring orbit, and knowing that the 
EOB formalism predicts the existence of a formal analog 
of the usual Schwarzschild light ring at the EOB dinien- 
sionless radius tlr = I/mlr, defined as the solution of 



i(uLR) =0, 



(6.55) 



with A{u) defined in Eq. (j5.2ip . it is natural to expect 
the (unknown) exact two-body version of the electric- 
quadrupole A potential to mathematically exhibit an 
analogous pole-like behavior of the form ~ (1 — tlr u)~^. 
As we shall discuss elsewhere, such a mathematical be- 
havior, linked to conside ring (within the EOB-simplifying 
approach advocated in [l8|) what would happen if one 
formally considered (unstable) circular orbits with u — >■ 
wlr, does not mean that there is a real physical singu- 
larity in the EOB dynamics near u = ulr, but it indi- 
cates that higher-than-2PN contributions to the electric- 

" (2) 

quadrupolc amplification factor electric (^) ~ 1 + u + 
ttj u'^+a'i^ + - • • will probably be slowly convergent, and 
will tend to further amplify the corresponding tidal inter- 
action. Such an extra amplification might, for instance, 
be physically important in the last orbits of comparable- 
mass neutron-star binaries (which will reach contact for 
values of u smaller than ulr)- 

This leads us to suggest that a more accurate value (for 
u < ulr) of the electric-quadrupole amplification factor 
is the following "resummed" version of Eq. 



i(2) 

1 electric 



''LRU 



3.56) 



where a° and are given by Eqs. (|6.9|) and (|6.10p . and 
where tlr =1/ wlr is the solution of Eq. (|6.55l) . Simi- 
lar resummed versions of the other amplification factors 
can be defined by incorporating in their PN-expanded 
versions the formal light-ring behaviors exhibited by the 
exact test-mass results (|6.44l) - (|6.47p . 

Let us finally discuss several possible approximate val- 
ues for in the proposed Eq. (I6.56|) . The simplest ap- 
proximation consists of using the "Schwarzschild" value 



^LR ~ ^- However, a better value might be obtained by 
taking a solution of Eq. (|6.55l) that incorporates more 
physical effects. This might require solving Eq. (|6.55p 
numerically, with A(u) being the full A potential (con- 
taining both Pade-resummed two-point-mass effects and 
the various tidal contributions). In order to have a feel- 
ing for the modification of tlr brought by incorporating 
these changes, let us consider solving Eq. (j6.55p when us- 
ing the following approximation to the full A potential: 
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Here the term + 2vu^ is the 2PN-accurate point-mass 
modification of Aiu), while the term —k is the leading- 
order tidal modification. Note that they have opposite 
signs. The corresponding expression of A{u) reads 



A, 



(6.59) 



The corresponding value of wlr =1/ ?^lr is the solution 
close to 1/3 of the equation 



MLR 



LR 



(6.60) 



If we could treat both 1/ and k as small deformation 
parameters, this would imply that, to first order in 
these two deformation parameters, the value of ulr(^', k) 
would be obtained by inserting the leading-order value 
mlr — 1/3 in the right-hand side of Eq. (|6.60p . This 
would yield 
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and 
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f-LR{v,K) = 3 



1-3^^ 



(6.62) 

Note that while comparable-mass effects (oc v) have the 
effect of decreasing tlr, tidal ones (oc k) have the op- 
posite effect of increasing tlr. Let us focus on the tidal 
effects, and consider the equal-mass case with Ri = R2 
One has a first order increase of tlr 



and k 
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equal to 
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where Ci = Grrii/ {c^ Ri) denotes the common compact- 
ness of the two neutron stars. This simple approximate 
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analytical formula shows that 5 f lr is very sensitive 
to the value of the compactness of the neutron star. If 
Ci = 1/6 = 0.166667, i.e. Ri — dGmi/c^ (roughly corre- 
sponding to a radius of 12 km for a 1.4 Mq neutron star), 

then <5"^^i tlr = 1.44 (/tpVO-09) wih be of order 1. [The 
value k^' = 0.09 being typical for Ci = 1/6; see, e.g., 
Table II in [|.] On the other hand, if Ri = Ri cV(Gtoi) 
is slightly smaller than 6, (5*"^^'rLR will quickly become 
much smaller than 1, while if i?i is slightly larger than 6, 

^tidal 

fLR will quickly become formally large (thereby in- 
validating the first-order analytical estimate (|6.63p which 
assumed Jtlr <C 3). These rough estimates indicate 
that, in many cases, tidal effects on tlr will be quite 
important and will significantly increase the numerical 
value of tlr. Note that an increased value of tlr will, 
in turn, increase the effect of the conjectured resummed 
2PN contribution a^u^/{l - tlru) to ifeioctric(")- 



VII. SUMMARY AND CONCLUSIONS 

Using an effective action technique, we have shown 
how to compute the additional terms in the reduced 
(Fokker) two-body Lagrangian L(yi, y2, 2/1, 2)2) that are 
linked to tidal interactions. Thanks to a general prop- 
erty of perturbed Fokker actions (explained at the end 
of Section 2, see Eq. (|2.20|) ). the additional tidal terms 
are correctly obtained (to first order in the tidal per- 
turbations) by replacing in the complete, unreduced ac- 
tion 2/1, 2/2] the gravitational field by the so- 
lution of Einstein's equations generated by two struc- 
tureless point masses mi, yi; 7712,2/2. This allowed us 
to compute in a rather straightforward manner the re- 
duced tidal action at the second post-Newtonian (2PN) 
fractional accuracy, by using the known, explicit form 
of the 2PN-accurate two-point-mass metric [36l - t39j . The 
main technical subtlety in this calculation is the regular- 
ization of the self-field effects associated with the com- 
putation of the various nonniinimal tidal-action terms 
~ / dT{Raf_,i3i, u'^)'^ + where, e.g., Raf^puix; 2/1,2/2) 
is to be evaluated on one of the worldlines that generate 
the metric (/^^ (so that i?Q^/3i/(2/i; 2/i, 2/2) is formally infi- 
nite). We explained in detail (in Section 3) one (algorith- 
mic) way to perform this regularization, using Hadamard 
regularization (which is equivalent to dimensional regu- 
larization at the 2PN level) . We then computed the regu- 
larized values, obtainable from the 2PN metric, of several 
of the geometrical invariants entering the nonminimal 
worldline tidal action terms. See Eqs. (I4.4l) - (|4.10p for 
the 2PN-accurate Lagrangians (for general orbits) of the 
three leading tidal terms (electric quadrupole, electric oc- 
tupole and magnetic quadrupole). We then focussed on 
the physically most useful information contained in these 
actions, namely the corresponding contributions to the 
effective-one-body (EOB) main radial potential, A{u), 
with u = G{mi+m2)/{c'^ rEOB)- Our Eqs. ((5?TO| . 
(|5.28p gave the explicit transformation between the pre- 



viously derived harmonic-coordinates tidal Lagrangians 
and their corresponding contributions to the EOB A po- 
tential. Using this transformation, we could finally ex- 
plicitly compute the most important tidal contributions 
to the EOB A potential to a higher accuracy than had 
been known before: namely, we computed the quadrupo- 
lar {i = 2) and octupolar {£ = 3) gravito-electric tidal 
contributions to 2PN fractional accuracy, i.e. with the 
inclusion of a relativistic distance-dependent factor of 
the type u'^^+'^{l + aiu + 02 m^), see Eqs. ([gl^ - lpTH)) 
and (|6.18P " (|6.22p . We also computed the quadrupolar 
gravito-magnetic tidal contribution, as well as a newly 
introduced time-differentiated electric quadrupolar tidal 
term, to IPN fractional accuracy, see Eqs. (|6.25p - (|6.28l) . 
(|6.30p - (|6.33p . Of most interest among these results is 
the obtention of the 2PN coefficient 012 entering the 
distance-dependence of the electric quadrupolar term. 
We found that this coefficient, Eq. (|6.10p . is always pos- 
itive and varies between H-3 and + 15.16071 as the mass 
fraction Xi — mi /{mi -I- 7712) of the considered tidally 
deformed body varies between and 1. In the equal- 
mass case, 7771 = ^772, i.e. Xi = ^, we found that 
= 6.07143. This value shows that, when the neu- 
tron star near their contact, 2PN effects are comparable 
to IPN ones. Indeed, contact occurs when the separa- 
tion r ~ Ri+ R2 = Gmi/{c^Ci) + Gm2/{c^C2) [where 
Ca = GmA/{c^ Ra), a ~ 1,2, are the two compact- 
nesses]. In the equal- mass case (withCi = C2), this shows 
that, at contact, u = G(mi+m2)/{c^ r) is approximately 
equal to itcontact — Ci. If we consider as typical neutron 
star a star of mass 1.4 and radius 12 km, we expect 
Ci ^ 1/6, i.e. Ttcontact ~ 1/6. The successive PN con- 
tributions to the distance-dependent amplification factor 

^i^cioc'tric(") = 1 + a" It -I- Q!2 of the electric quadrupo- 
lar tidal interaction for the first body become then, at 
contact. 
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(7.1) 



where one sees that the 2PN (0(u^)) contribution is nu- 
merically comparable to the IPN one. This suggests that 
the PN-expanded form of the tidal amplification factor 
„{u) is slowly converging and could get compara- 



(2) 

1 clcctr 



A 

ble, or even larger, contributions from higher powers of 
77 (i.e. 3PN and higher terms). In order to get a feeling 
about the possible origin of this slow convergence of the 
PN expansion, we followed the approach of 5J|, i.e. we 
looked for the existence of a nearby pole (in the complex 
77 plane) within the formal analytic continuation of the 

considered function 4i^oioctric('")- I^^^- 01 considered 
the energy fiux F as a function oi x = (GM fi/c'^)^/'^; 
it pointed out that F{x) had (in the test-mass limit) a 
pole at the light-ring value x = 1/3, and recommended 
to improve the PN expansion of F{x) (for x < 1/3) by 
a Pade-type resummation incorporating the existence of 
this pole in F{x).] By computing the exact test-mass 
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^ (2) 

limit of the function electric ("")' found that it did 
formaUy exhibit a pole located at the light-ring value 
ytestmass ^ 2^/3^ (lOI)) . Such a pole is also present 

in other amplification factors, see Eqs. (|6.45p - (|6.47l) . and 
we discussed its origin. [Note that two equal-mass neu- 
tron stars will get in contact before reaching this pole. 
However the idea here is that the "hidden" presence of 
this pole in the analytical continuation of the function 

^ieiectric(") behind the bad convergence of the Tay- 
lor expansion of this function in powers of u.] This led 
us to suggest that one might get an improved value of 
the tidal amplification factors electric formally 
incorporating the presence of this pole in the following 
"Pade-resummed" manner 

2 

^llctric(") = 1 + + «^ I , (7.2) 

where tlr = 1/ulr is the (EOB-defined) dimensionless 
light ring radius, i.e. the solution of Eq. (|6.55p . with A{u) 
defined by Eq. Let us point out that Eq. ((7?^ is 

equivalent to saying that the 2PN coefficient a2 becomes 
replaced by the effective distance-dependent coefficient 
af{u) = - TLRu). Note that af{u) > a^. In 

particular, for the "typical" compactness Ci = C2 ^ 1/6 
considered above, and when using the unperturbed value 
of tlr, i.e. = 3, the effective value a2^{u) will, at 
contact (i.e. when u — Wcontact — Ci ^ 1/6), be equal to 
af (ucontact) ^ ai/(l-3Ci) - ai/(l-3/6) ^ 2ai ~ 12. 
We recalled in the Introduction that several comparisons 
between the analytical (EOB) description of tidal effects 
and numerical simulations of tidally interacting binary 

neutron stars [5i, i9i, ilO| have suggested the need for sig- 

^ (2) 

nificant amplification factors ^1 eiectric(")' parametrized 
by rather large value of ctj. However, up to now, the nu- 
merical results that have been used have been affected by 
numerical errors that have not been fully controlled. In 
particular, in the recent comparisons P, [l3|j one did not 
have in hand sufficiently many simulations with different 
resolutions for being able to compute and subtract the 
finite-resolution error. We hope that such a more com- 
plete analysis will be performed soon. We recommend to 
compare resolution-extrapolated numerical data to the 
pole-improved amplification factor (|7.2p . As discussed 
in Section IVIl it might be necessary to use as value of 
^LR the improved estimate obtained from the full (tidally 
modified) value of the A potential. This suggests (espe- 
cially for compactnesses Ci < 1/6) as discussed above 
that tlr might be significantly larger than 3, thereby 
further amplifying the effective value of ctg during the 
last stages of the inspiral. 

The present study has focussed on the 2PN tidal ef- 
fects in the interaction Hamiltonian. There is also a 2PN 
tidal effect in the radiation reaction which has contri- 



butions from various tidally-modified multipolar wave- 
forms. The tidal contribution to each (circular) multipo- 
lar gravitational waveform can be parametrized (follow- 
ing Refs. @,I3) as an additional term of the form 

hf:L^\:>:) - E '^^i' ^^i' ™(-) ' (7-3) 
J 

where: x = (G(mi -I- m2) il/c^)^/'^, J labels the vari- 
ous tidal geometrical invariants, such as Jq = Gq,^ G"^ , 

hg^^'^{x) denotes the leading-order (i.e. Newtonian- 
order) tidal waveform, h''i^^^^^{x) the effect of tails 
[55I [Sa j and their resummed EOB-form [s^l ; while 

/,(;2™(x) = 1 + + pi'''-^x' + ... (7.4) 

denotes the effect of higher PN contributions. The IPN 
coefficient /Sj-^"^^^ is known 0,[ll. The other IPN coeffi- 
cients needed for deriving a 2PN-accurate fiux can be ob- 
tained from applying the simple IPN-accurate formalism 

of [40]. It is more challenging to compute the 2PN coef- 
( J 22) 

ficient . Indeed, this requires applying the 2PN- 

accurate version [56j of the Blanchet-Damour-Iyer wave- 
generation formalism [58l - l60| to the tidal-modified 
Einstein equations ()2.13p . Let us, however, note that, 

though, from a PN point of view, the 2PN coefficient 
(J 22) 

1^2 ° contributes to the phasing of coalescing binaries 
at the same formal level as the dynamical 2PN coefficient 
Q!2 determined above, it has been found in [l^ that 

(if /3^-^"^^^ ~ ai) it has a significantly smaller observable 
effect. 

Let us finally point out that our general result (|2.20p 
opens also the possibility to compute the 3PN coefficient 

ftg in the PN-expanded amplification factor of the electric 

^ (2^ 

quadrupolar tidal interaction: ^1 eiectric(^) ^ 1 + ct? u + 
a2 -\-0{u'^). This computation would, however, 

be much more involved than the calculation of 02 because 
of the technical subtleties in the regularization of self- field 
effects at the 3PN level [i^, l6ll-{63l that necessitate to use 
dimensional regularization [2^, [2g] instead of Hadamard 
regularization. 
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